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The Shape of the Strongest Column 


JOSEPH B. KELLER 


Communicated by C. TRUESDELL 


Abstract 


The problem of determining that shape of column which has the largest 
critical buckling load is solved, assuming that the length and volume are given 
and that each cross section is convex. The strongest column has an equilateral 
triangle as cross section, and it is tapered along its length, being thickest in 
the middle and thinnest at its ends. Its buckling load is 61.2% larger than that 
of a circular cylinder. For columns all of whose cross sections are similar and 
of prescribed shape—not necessarily convex—the best tapering is found to 
increase the buckling load by one third over that of a uniform column. This 
result, which was independently obtained by H. F. WEINBERGER, is originally 
due to CLAUSEN (1851). For a uniform column, triangularizing is shown to 
increase the buckling load by 20.9% over that of a circular cylinder. The results 
lead to isoperimetric inequalities for the buckling loads of arbitrary columns. 


1. Introduction 

Of all columns having prescribed length and volume and convex cross sections, 
which is the strongest, 7.e., which has the largest critical buckling load? We 
shall solve this problem and show that the strongest column is not uniform but 
is tapered, being thickest at its center and thinnest at its ends. We shall also 
show that the cross section of the strongest column is not a circle but is instead 
an equilateral triangle. Appropriately tapering a column of circular cross section 
increases its critical buckling load by one third over that of a uniform column. 
Changing the cross section from a circle to an equilateral triangle increases the 
critical buckling load by 20.9%. Thus the combined modifications of tapering 
and triangularizing increase the critical buckling load by 61.2% over that of a 
circular cylinder. 

Our principal result can also be stated as an isoperimetric inequality. In 
this form it states that the critical buckling load of any column is less than a 
given function of the length, volume and Younc’s modulus of the column. The 
given function is just the critical buckling load of the strongest column with 
the same length, volume and Younc’s modulus. In the course of our analysis 
two other isoperimetric inequalities arise. One pertains to an untwisted tapered 
column with similar cross sections and the other to cylindrical columns. The 
former can also be stated as an isoperimetric inequality for a certain eigenvalue 
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problem, without reference to columns. The latter is related to a geometrical 
problem which is not completely solved. In the last section we present some 
results for the nonlinear buckling problem. 

For the sake of clarity we first determine the strongest cylindrical column 
of arbitrary convex cross section. Then we determine the strongest untwisted 
column, all cross sections of which are similar. Upon combining these two results 
we find the strongest untwisted column with convex cross sections. Finally we 
show that this column is also strongest among all columns with convex cross 
sections. This procedure is also economical since we use the simple results in 
the final section. 

Our analysis is based upon the six equations of equilibrium of a column of 
arbitrary shape. To these we add expressions for the three stress couples in 
terms of the curvatures and twist of the column. These three relations are the 
usual generalizations of the Bernoulli-Euler theory. We then seek solutions of 
the nine equations subject to the conditions that the stress couples vanish at 
the ends of the column and that the applied load have a prescribed value. To 
determine the critical buckling load we linearize the equations about the straight 
unbuckled state. The linearized problem is an eigenvalue problem and has non- 
zero solutions only for special values of the load. The smallest of these is the 
critical buckling load. It depends upon the shape of the column. We next apply 
variational techniques to determine that shape which maximizes the critical 
buckling load. In this way we obtain the results described above. 

I wish to thank CLIFFORD TRUESDELL for proposing this problem to me. 
He also suggested it to HANS F. WEINBERGER, who independently obtained the 
results of Section 3 on untwisted columns with circular cross sections. However, 
as TRUESDELL later pointed out, that result had been obtained by CLAUSEN 1 
in 1851. LAGRANGE? in 1773 first treated that problem but arrived at the wrong 
result due to computational errors?. 


2. Cylindrical columns of arbitrary convex cross section 


For simplicity let us begin by considering cylindrical columns of arbitrary 
cross section. Let L and V denote the length and volume of such a column and 
A=V/L its cross sectional area. Let the x-axis pass through the centroid of 
each cross section of the column in its straight or unbuckled state. Suppose 
that I’ is the largest moment of inertia of the cross section about a line in its 


' Uber die Form architektonischer Saulen. Bulletin physico-mathématique de 
YAcadémie 9, 368—379 (St. Pétersburg 1851). Also Mélanges Mathématiques et 
Astronomiques, Tome I (1849—1853), pp. 279—294. St. Petersburg 1853. This work 
1s summarized in TopHuNTER, I., & K. Pearson, A History of the Theory of Elasticity 
and of the Strength of Materials, vol. 2, pp. 325—329. Cambridge 1893. 

» Sur la figure des colonnes. Miscellanea Taurinensia (Royal Society of Turin) 
Toames VW, (IAS, MAO—N77B, AGO Oeuvres, vol. 2, pp. 125—170. Gauthier-Villars 
1868. Summarized in TopHUNTER & PEARSON, vol. 1, pp. 66—67. Cambridge 1886. 
A more complete early history of this and related problems is given by TRUESDELL 
The Rational Mechanics of Flexible or Elastic Bodies, 1638—1788, L. Euleri Opera 
Omnia II 11, (in press). 


aris ‘ pert : 
; These errors are pointed out in his Oeuvres by the editor. They are not mentioned 
in TODHUNTER & PEARSON. 
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plane through the centroid and that the line is the y-axis. Similarly let the 
minimum moment of inertia be J. Then the corresponding line is the z-axis. 

Let us consider a state of equilibrium of the column in which a compressive 
load Yj is applied longitudinally at its ends. Let us further suppose that the 
deflection of the column from the straight unbuckled state is small. Then the 
equilibrium is governed by linear equations and every state is a superposition 
of a deflection in the xy-plane and another in the xz-plane. These deflections 
y (x) and z(x) satisfy the equations 


Ly 
Veter ag 2) (Oli VAL) = 0, (1) 


z=0, 2(0)=2(L)=0. (2) 


The boundary conditions in (1) and (2) apply if the column is pinned at its ends, 
and £ denotes the YounG’s modulus of the column material. 

The least positive value of 7) for which (1) has a non-trivial solution is 
mE IL while for (2) it is m?E I'L. The critical buckling load is the minimum 
of these two values. Since by assumption J<J/’, we have for the critical load 


Ty = 2 ETL. (3) 


The strongest column is that for which 7) is as large as possible. Since E, L 
and V, and therefore A, are prescribed, the strongest column is one for which I 
is largest. But J can be made arbitrarily large if the cross section is nonconvex. 
Therefore the problem of finding a strongest column makes sense only if the 
cross section is required to be convex. Then the cross section of the strongest 
column is the solution of the following geometrical problem. Problem: Of all 
plane convex domains of area A, which has the largest value of J? Here I 
denotes the minimum value of the moment of inertia of the domain about any 
line through its centroid. 

This problem does not seem to have been investigated previously. It is 
easy to show that the domain which solves the problem has the same moment 
of inertia about every line through its centroid which lies in the plane. However 
this domain is not a circle. This follows from the facts that the circle has the 
smallest polar moment of inertia of all plane figures of given area and that I 
is one half the polar moment for any figure having the same moment J about 
every line through the centroid. Every regular polygon has the same moment 
about every line, through the centroid, and of these the equilateral triangle has 
the largest J for given area. Therefore it has been conjectured by PETER UNGAR 
that the equilateral triangle is the solution of the above problem. GEORGE STELL 
has proved that it is the solution among all figures having an -fold symmetry 
axis with »=3. The complete solution has not yet been given, but there is 
little doubt that it is the equilateral triangle. If this is the case, we have the 
isoperimetric inequality for convex domains 


ee A2(096.225)i (4) 
6/3 


Equality holds for the equilateral triangle. 
AO 
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By using (4), (3) yields the isoperimetric inequality for the critical buckling 
load of cylindrical columns with convex cross sections 

2B A EA? 

SEE ete 5 

VSR Es (.949 703) (5) 

Equality holds for the equilateral triangular cross section. The numerical co- 

efficient in (4) is 20.9% larger than the value 7/4=.785 for the circular cross 

section. 
3. Tapered untwisted columns with similar cross sections 


Let us now consider tapered columns. We shall still assume that all cross 
sections are similar and that the principal axes of inertia are parallel. Thus 
the column is not twisted. However the area A(x) of the cross section may vary 
with x subject to the condition 


1 
f A(x) dx=V. (1) 
0 
The moments J and /’ are related to A by the equations 
LG) = ANA yk Ate) (2) 


The constants « and «’, with «<«’, depend upon the shape of the cross section. 
Now (1) becomes 


Yost peat Y= ¥(0)=¥(L) =O. (3) 


We seek that function A(x) satisfying (1), which maximizes the lowest eigen- 
value 7) of (3). This function will be the area variation of the strongest column. 

Before solving this problem we first introduce the new variable = x/L and 
the eigenvalue parameter A= 7) L?/x E and consider y and A as functions of &. 
Then (3) and (1) become 


y’+AA(E\y=0, y(0)=y(1)=0, (4) 
1 


J At) CEN Ts (5) 


Now we assume that A(&) is the solution. We introduce a family of functions 
A(&, e) depending differentiably on a parameter e, satisfying (5), and containing 
the solution A(&)=A(E,0). Then for each e¢, (4) will have a solution y(6é, e) 
with the corresponding lowest eigenvalue A (e), both of which depend differentiably 
on ¢. Next we differentiate (4) and (5) with respect to ¢, obtaining 


Ye tAA*y,=21A,ASy—f,A%y, (0) =¥,(1) =0, (6) 
1 
JA. de =0; (7) 


The derivative y, satisfies an inhomogeneous linear differential problem (6) 
and y satisfies the corresponding homogeneous problem (4). Therefore the right 
side of the differential equation in (6) must be orthogonal to y. At e=0, A,=0 
since by assumption A is a maximum at e=0. The orthogonality condition is 
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then 1 
ASA tan) d&=0.., (8) 


Since (8) must hold for every A,(&) satisfying (7), it follows that A-%y? is a 
constant which we denote by c!. Thus for the strongest column we have 


yr = cb AB, (9) 
To find y we first eliminate A from (4) by means of (9) and obtain 
y’+Acy*t=0, y(0)=y(1)=0. (10) 
The solution of (10) is 
pS Ee Ay ceyh e 
Here yy is given by yh — a (Ac)A. (12) 


Upon replacing y in (11) by A, with the aid of (9), we obtain the following equation 
for the solution A(&): . 


é =|= ; * = (4) }- sin Cae (13) 


The constant A4)=A(#) in (13) can be determined as a function of V/L by in- 
serting (13) into (5). This yields 


3 Ay 
VIL =o PAeyde a2 fedd = 340, 1 
ai (é) a : (14) 


To find A let us use (9) to eliminate y from (4), which becomes 
AA" +4(A’)??+21=0. (15) 


From (13) we find that A’($)=0 and A’’($)=2?A,/2. Upon inserting these 
values into (15) and using (14) for Ay we obtain 
47? (V\2 
Nerd ae ee, 
Now from the definition of A we find the critical buckling load of the strongest 
column. The result may be stated as an isoperimetric inequality 


Ty *F a VOLE. (17) 


Equality holds only for the strongest column of given volume, the shape of 
which is given by (13). For it the critical buckling load is $ times as large as that 
of a cylindrical column of the same volume, length and cross section. 

Let us now combine the results of this and the preceding section. Then we 
find that the strongest untwisted column, all of whose cross sections are similar 
and convex, is tapered according to (13) and has an equilateral triangular cross 
section. Upon using the value «= (61/3) appropriate to this triangle, we obtain 
from (17) the isoperimetric inequality 


The ce ee a ae (1.26627). (18) 
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Equality holds only for the column just described. For it the critical buckling 
load is 1.612 times as large as for a circular cylindrical column. In the remaining 
sections we shall show that (18) also holds for twisted columns, without the 
restriction that the cross sections be similar. 

Our result (16) can be expressed as an isoperimetric inequality for the eigen- 
value problem (4) without reference to columns. It can be formulated thus: 
Let A be the lowest eigenvalue of the problem 


y+ Ag(x)y=0, y(0)=y(L)=0, 9{%)>0. (19) 
Then ee S 
Assia | fate) dx. (20) 
Equality holds if q(x) =A?(«/L), where A is given by (13). 
3 
7) 


me 
—_— _-—_——= 


T I | | 
0 B7 DR UF to Oban S aN OB Pal ae Tee ee 10 


Fig. 1. The area (solid curve) and square root of the area (dashed curve) of the strongest column as functions of distance 


1 
along the column. These curves are based upon equations (21) and (22). The vertical coordinate is7=A L/V or n= (AL/V)?, 
and the horizontal coordinate is =%/L. Here A(x) is the cross sectional area at the distance x along the column, L is 

the length and V the volume of the column 


The equation (13) for the shape of the strongest column can be simply ex- 
pressed in terms of a parameter # which varies from 0 to x. Then, as CLIFFORD 
TRUESDELL has observed, (13) and (14) yield 


me es ae 
Aves cin sin? ?, (21) 
1 Wee 


4. Arbitrary columns with convex cross sections 


Let us now consider any column of length L and volume J, all cross sections 
of which are convex. Among them we seek the strongest one. To find it we 
begin by considering any such column in a state of equilibrium in the absence 
of external distributed forces or couples. The vanishing of the net force and 
torque on each portion of the column leads to six equations [LovE p. 387 eq. 10, 


Papesieqs 44 \+ A SNe eo 1) 
N, —Tx+Nr=0, (2) 
Re NN eet (3) 
G,—G't+Hx'—N'=0, (4) 
G,— Hx+Gr+N=o0, (5) 
H,—Gzx'+G'x=0. (6) 


1 Love, A. E.H.: A Treatise on the Mathematical Theory of Elasticity, 4t ed. 
Cambridge Univ. Press. 1927. 
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In these equations N and N’ are the shearing forces, T is the tension, G and G’ 
are the flexural couples, H is the torsional couple, x and x’ are the components 
of curvature of the central line and 1 is the twist of the rod. The subscript s 
denotes differentiation with respect to arclength along the central line. The 
forces and couples are expressed as components along the principal torsion- 
flexure axes of the rod. 

In the usual approximate theory of a naturally straight thin rod, which is 
a generalization of the Bernoulli-Euler theory, it is assumed that the stress 
couples are related to the curvatures and twist of the rod by the equations 
[LovE p. 388, eq. 12, p. 397, eq. 28]! 


G=EIx, (7) 
G=ETI’x', (8) 
H =C(t — %). (9) 


In these equations E is the YounG’s modulus of the rod material, which is 
assumed to be isotropic, J and J’ are the moments of inertia of the cross section 
of the rod about its principal axes, C is the torsional rigidity and T, is the twist 
in the unstressed state. The nine equations (1)—(9), together with appropriate 
boundary conditions, suffice to determine the nine quantities N, N’, T, G, G’, 
Til be. ad Bane), TH 

We wish to consider the buckling of a naturally straight column subjected 
to a compressive axial load. In the straight or unbuckled state x=x’=0. Then 
from (7) and (8) it follows that G=G’=0. But then (4) and (5) yield N= N’=0. 
Therefore (3) and (4) become 7; =H,;=0 or T= — 7)=constant, H=constant. 
We shall assume that no twisting couple is applied at the ends of the rod. Con- 
sequently H =O, and then (9) shows that t=7,). Thus in the straight state only 
T and t may differ from zero. If we consider a buckled state which is sufficiently 
close to the unbuckled state, we may assume that all quantities except T and t 
are small while T and 7 are close to their values in the unbuckled state. Then 
we may linearize (1) —(9), which become 


N;— N’1,— 1g’ =0,; (10) 
Ne + To% tN % =0, (11) 
L,=0, (12) 
G.— Gita N= 0, (13) 
G+ Ga N=0, (14) 
fia) (15) 


From (15) and the boundary condition that no couples act at the ends, it follows 
that H=0, and (9) still shows that t=t,). Similarly (12) yields 7 = constant, 
and if J) is the applied load, then T= — 7). 

The linearized buckling problem is that of solving (10), (11), (43), (14), (7) 
and (8) for N, N’,G,G’,x and x’ when I, I’ and ty, are given functions of s, 
and E is a given constant. At the endpoints s=0 and s=L we require that 


1 See the previous footnote. 
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no couples act. Thus Go) —6'(0) = G(L) = GD =0. (16) 
This is an eigenvalue problem and will have non trivial solutions only for certain 
special values of the constant 7). The smallest buckling load is the least value 
of Ty for which a non trivial solution exists. 


If we eliminate x and x’ from (10) and (11) by means of (7) and (8), we obtain 


N,—N'%p— 8 G20, (17) 
Ni 4+ N1+2G=0. (18) 


We must now solve (13), (14), (17) and (18) for N, N’, G and G’ subject to (16), 
which is possible only if Jj is an eigenvalue. We now seek that column, charac- 
terized by t,, J and I’, for which the first eigenvalue is largest. It must have 


the volume V, so L 
f A(s)ds=V. (19) 
0 


5. The variational problem 


Let us now consider a one parameter family of columns depending differen- 
tiably upon the parameter ¢ and containing the strongest column at e=0. Thus 
T, J and I’ are all differentiable functions of ¢ as well as of s. Then the solution 
of the buckling problem will also depend differentiably on e. Then we obtain 
from (13), (14), (17), (18), (16) and (19) of Section 4 respectively, with 6/de 
denoted by a dot, the variational equations 


Gt, NSO) (1) 
Git Gri == Gt; (2) 
N4N1— $2. G=—Nvi—(2)6, (3) 
N— Way + 28,6 =N'in + (F276, (4) 
G (0) =G'(0) =G(L) =G(L) =0, (5) 

faas=o. (6) 


The equations (1) (5) form an inhomogeneous system for the four derivatives 
G,G’, N and N’. The corresponding homogeneous system is just that satisfied 
by.G,.G, N, and N’,. i.¢., (43); (44), (47), (18) and (16) of Section 4. Therefore 
in order that the inhomogeneous system have a solution, the right sides of (4) 
to (4) must be orthogonal to all the solutions of the adjoint homogeneous system. 
The only solution of this system is (—N’', N, —G’, G). Thus the orthogonality 
condition is 

L 


f[-2W@+NG)%—6? (3-)-@ (e)| ds=0. (7) 
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In (7) we choose I’=I=0 and assume that Tr=0 at e=0. Then since (7) must 
hold for every choice of t), we conclude that for e=0 


ZG NG (8) 


Let us now deduce some conclusions from the foregoing equations. First we 
multiply (13) of Section 4 by —G’ and (14) of Section 4 by G, add and use (8) 
to obtain 

G,G — G’G, + (G24 G) m=0. (9) 


We rewrite (9) and define f, as follows: 


(GIG). = 
(G'/G)?-+1 a Co Salis (10) 


Upon integrating (10) and adding an appropriate constant to f, we have 


G’'/G = — tan. (14) 
Let us define g by 
ge A Eels AEE (12) 
Then from (41) and (12) we find 
G == gsinf, ..G=g.cosp: (13) 
Now we define 1 by 
n? = N’2 4+ N2, (14) 
Then from (8) and (14) we get 
N2=ncosp, N=nsinp. (15) 
Now we multiply (13) of Section 4 by G, (14) of Section 4 by G’ and add, obtaining 
3 (G?+ G’*),—GN’'+G’'N=0. (16) 
By using (13) and (15) in (16) we find 
Bee te O. (17) 
Hence, if g=-0, we have 
f= 0. (18) 
We now insert (13), (45) and (18) into (17) and (18) of Section 4 and obtain 
1) 1 = 
(é..+ 97 e)esinp=0, (19) 
+ 3,8] ¢005 B= 0. (20) 
8ss ETI 
If g sin B 0 and g cos B 0, we conclude from (19) and (20) that 
eau (21) 
ee (22) 


From (21) and (2) of Section 2, with «=«(s), %’=«'(s), it follows that a=’, 
If we set A= Ta E, we may rewrite (22) in the form 


85s tAA*g=0. (23) 


Pe 
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From (16) of Section 4 and the definition of g, we have 
gO) = 6 (ENO; (24) 


If 2 denotes the lowest eigenvalue of the problem (23), (24), then Ty>=aEFA. 
Thus for a fixed A®, J, is maximized by maximizing «. Since <=’, the problem 
of maximizing « is just the same geometrical problem considered in Section 2. 
Then since the best value of « is a constant, independent of s, the problem of 
finding A(s) to maximize 4 is just that solved in Section 3. Thus if sin B=:0 
and cos B=£0, the strongest column is that described toward the end of Section 3. 

The exceptional case in which sinf=O or cosf=0 occurs if ty=0, 1.e., if 
the column is untwisted. Then we must compare the lowest eigenvalues for the 
two cases B=0 and z/2. The smaller of the two is the critical buckling load. 
The former case yields (22)—(24), and the latter yields the same equation with 
a’ in place of «. Thus the critical buckling load is the smaller of «EA and «’ EA. 
The determination of the best cross section is thus just that of solving the geo- 
metrical problem of Section 2, and then the determination of the best A(s) is 
the problem of Section 3. Thus the strongest column is that found in Section 3. 


6. The nonlinear buckling problem 


In the preceding sections we considered the linear buckling problem. That 
problem was obtained by linearizing the original nonlinear problem around the 
straight unbuckled state. It yields the surprising result that there is no buckled 
state unless the load is equal to an eigenvalue. But we expect buckled states 
to exist for any load exceeding the lowest eigenvalue. In fact, if we consider 
deflections in a single plane, we expect ” buckled states to exist when the load 
exceeds the n" eigenvalue of the linear problem, but does not exceed the (m+ 1)*. 
These buckled states correspond to the m buckling modes associated with the n 
eigenvalues less than the load. This expectation is borne out for cylindrical 
columns for which the explicit solutions of the nonlinear problem (the elastica 
problem) are known. However for tapered columns, such as those considered 
here, it has not yet been shown that this is the case. The corresponding statement 
has been proved to be true by I. I. KoLopNER! for a somewhat similar problem. 

As a first step toward verifying this conjecture, we shall show that there 
is a buckled state for any load slightly larger than any buckling load. To do 
so we shall consider deflections of an untwisted column in a single plane and 
let ¥(s) denote the angle between the column and a fixed axis in this plane. 
Then, as follows from the equations (1)—(9) and (16) of Section 4, ®(s) satisfies 
the equations 


([3,),+Asin®=0, 8,(0) =%,(L) =0. (1) 


The constant A= T,E and the positive function I(s) are given, and @(s) is to 
be found. 


It is convenient to introduce the constant a and the function w defined by 


O=au, a=8(0). (2) 


1Comm. Pure Appl. Math. 8, 395—408 (1955). 
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Then w satisfies the equations 
(Tj Aa sin(au)— 0, ~u(0)—= 1 :74,(0)=0, 4, (B=0. (3) 
If we let a tend to zero in (3), we obtain the linear problem 
Ga) A¢=30,, UOl=—=4, “(OJ —0, 4(L)=0. (4) 


This problem has a discrete set of increasing positive eigenvalues J) and eigen- 
functions w?. 

Let us suppose that a is given and that w and / are both to be found. Then 
we assume that for a small 


U = Uy + OP ty boo, (5) 


a7) Aa? EE ey, (6) 


We have expanded in powers of a? because (3) is even in a. Insertion of (5) and 
(6) into (3) and a simple calculation yield the result 


L 
f (un)* ds 
an 
Mier (7) 
J (un)? ds 
0 


Since 4), >0, we can solve (6) for a in terms of 4 provided A= A?. If we denote 
the solution by a, (A), we have 


0 

a, (8) = (Set) + Ol — m)]. 8) 
By using this result, it is easy to prove that (3), and therefore (1), has a solution 
for A—A°* sufficiently small and non-negative. This shows that the nonlinear 
problem has buckled solutions for any load in a certain half-neighborhood of 
each buckling load of the linear problem. It also shows that the lowest buckling 
load of the linear problem is also the lowest buckling load of the nonlinear problem, 
at least for slightly buckled states. Of course the conjecture concerning the 
number of buckled states remains to be proved. 


The research reported in this paper has been sponsored by the Office of Naval 
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An Optimal Poincaré Inequality for Convex Domains 


L. E. PAYNE & H. F. WEINBERGER 


Communicated by C. TRUESDELL 


1. Introduction 


Let G be a convex u-dimensional domain with boundary C. It is easily seen 
that the lowest eigenvalue of the free membrane problem 


Av+tpwv=0 inG 


Ale 
EY dvJen=O0 onC 


is zero, the eigenfunction being any constant. 
This corresponds to the fact that the solution of the interior Neumann problem 


Ag=0: inG 
(1.2) | 
dgl/on  givenonC 


is only determined to within a constant. The latter is to be fixed by a normali- 
zation such as 


(1.3) le ae 


The authors have previously introduced a method for bounding the point- 
wise value as well as the Dirichlet integral of a solution @ of the exterior Neumann 
problem in terms of a boundary integral of (0/0)? [3]. In order to extend 
this method to the interior Neumann problem one needs a lower bound for the 
second eigenvalue mu, of (1). This eigenvalue is characterized by the minimum 
principle 
J | grad u|? dG 
G 


ie lb = mi 
ae) Ma fudge-o fwdG 
G G 
A lower bound for j,. can be used in the interior Neumann problem in the follow- 


ing manner (cf. [3]). Let f be a vector field which is piecewise continuously 
differentiable throughout G and points outward on C, so that 


(1.5) Push O ONG 


For example, if G is star-shaped with respect to the origin, we may take } to 
be the radius vector. By the divergence theorem and the inequality a2+b2=2ab 
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we have, if m is normalized by (1.3), 


= 


$e aaC= f [praivi+ 2pf- grad y] dG 
=! [div /+|f[*]@G + f|grad g|24G 
(1.6) 6 
< [1+ yg! max (div / + |f]?)] f | grad 9/246 
G 
= [1+ wat max (div f+ |f|?)] $y a/andC. 
(G 


Consequently by SCHWARz’s inequality 


> 


(1.7) Deiat C Ss.1 + yz’ max (div f+] /]2) \Ps n) (dg/on)2dC. 


If J° is a fundamental solution of LAapLace’s equation with its singularity 
at the interior point P, we have 


|p 
(1.8) < ee 7) (2g|on) acl Cae 


(F. 


$ (I’dg/on — y OLén) ) ac| 


+ [4+ w3*max (div /+|/|? {oe n) 1 (arjenj2acl’). 
Thus, if a lower bound for “7, is known, gy(P) may be explicitly bounded in terms 
of a the square integral of dg/dn. 

These results can be extended to general second order differential equations 
(cf. [3]). 

In this paper we shall show that for a convex domain G in any number of 
dimensions 
(1.9) a= ge DH 
where D is the diameter of G. This is the best bound that can be given in terms 
of the diameter.alone in the sense that w,D? tends to z* for a parallepiped all 
but one of whose dimensions shrink to zero. 

The inequality (1.9) is in general false for non-convex domains. In fact, 
for a sequence of domains which tends to two disjoint subdomains, mw, tends 
to zero. For the special class of domains G which are symmetric about all the 
coordinate planes of a rectangular coordinate system and have the property 
that the intersection of G with any line in a coordinate direction is simply con- 
nected, the authors have previously obtained an inequality of the form (1.9) 
with D replaced by the maximum length of intersection of G with a line in any 
coordinate direction [4]. 

A simple upper bound for jw. for any m-dimensional domain G in terms of its 
volume V is given by the isoperimetric inequality 


he od 
(1.10) Peay, bela”, 


288 L. E. Payne & H. F. WEINBERGER: 


where K, is the volume of the unit 2-sphere and #,, is the lowest positive root 
of the equation 


(1.44) b Iyn(P) — (2% — 1) yn (b) = 0. 


Equality is attained when G is a sphere. 

For n=2 this inequality was conjectured by KORNHAUSER & STAKGOLD [2] 
and proved by SzEG6 [5]. For general » the proof was given by one of the 
authors [6]. 

The eigenvalue j itself is of interest in a variety of problems arising in 
mathematical physics. In two dimensions it is proportional to the square of 
the cutoff frequency of the lowest H-mode of a wave guide [2]. In three dimensions 
it is proportional to the lowest resonant frequency of an acoustic resonator with 
perfectly rigid walls. It is also inversely proportional to the relaxation time for 
diffusion in a body with perfectly reflecting boundary. 

The proof of the lower bound (1.9) is based upon a lemma concerning a class 
of Sturm-Liouville systems. This lemma, which is of some interest in itself, is 
stated and proved in § 2. The inequality (1.9) is proved for two dimensions in 
§3 and for higher dimensions in §4. 


2. A one-dimensional lemma 
In order to prove the lower bound (1.9) we require a somewhat stronger 
version of its one-dimensional analogue. It is the following lemma. 
Lemma. Let p(y) be a non-negative convex function of y defined on the interval 


OSySL; then for any piecewise continuously differentiable function u(y) which 
satisfies 


(2.1) AE ACI EON LS At 
at follows that 
(2.2) J ety) [u(y Pay = EN) P(y) (u(y) Pay. 


Proof. We assume for the moment that # is strictly positive and twice differ- 
entiable. Then the function v which minimizes the quotient 


16 L 
(2.3) [pw ?dy/fpurdy 
0 0 


among functions w satisfying (2.1) must satisfy the Sturm-Liouville system 
LZ, p. 348] 


0.4 [pu sr ap ves0; 


v'(0) =v'(L)=0, 


where 4 is the minimum value of the quotient (2.3). We divide the equation 
(2.4) by p, differentiate with respect to y, and introduce the new variable 


(2.5) w=v' pi, 


Poincaré Inequality for Convex Domains 289 


The function w satisfies the Sturm-Liouville system 


Ae 4 p”’ 3 p? 
w+ wtiAw=O, 
(2.6) Se og Ane 


w(0) =w(L)=0. 


Because of the convexity of p the term in square brackets is non-positive. Hence, 
multiplying (2.6) by w and integrating by parts, we obtain 


(2:7) | es 


Since w(0)=w(L)=0 the quotient on the right of (2.7) is bounded below by 
the first eigenvalue of the vibrating string with fixed ends. Thus 

(2.8) Ve at, sad Deas 

Since A is the minimum of the quotient (2.3), (2.2) is proved when # is strictly 
positive and twice differentiable. 


If % is any function defined on the interval OX y<L, the function 


L Ethie 
(2.9) uo) = ay) —|fody] Spaay 
will satisfy (2.1). Hence (2.2) implies 
L L L fel L 2 
(240)  fpa'*dyzaPL*) fpitdy—|fpay] | fpaay] |. 
0 0 0 0 


Clearly (2.10) is valid for the uniform limit of admissible functions #. In parti- 
cular, then, p may be any non-negative convex function of y. Thus the lemma 
is proved. 


Remarks. 1. The convexity of # was used only to show that the square bracket 
in (2.6) is non-positive. For this purpose it is sufficient to assume that p~? is 
a concave function of y. Therefore the lemma actually holds under this weaker 


condition. 


2. By the minimax theorem [J/, p. 352] we can show that if p~2 is a concave 
function of y, the eigenvalues of the Sturm-Liouville system (2.4) satisfy the 
inequality 
ei) Ay 2d )2 RL REMAND). ba 


3. Equality in (2.11) is obtained if and only if p~? is linear in y. If pis assumed 
convex, it must then be constant. 


3. The two-dimensional case 


Let G be a convex plane domain with boundary C. Let uw, be defined as the 
infimum! of the quotient 


(3.4) f|gradu|?dG/ [wdG 
G G 


1Tf the boundary C is smooth so that the problem (1.1) possesses eigenvalues, 
[ly is the second eigenvalue of (1.1). 
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among functions which have bounded second derivatives in G and satisfy 


(3.2) Uae ah 


Let « be such a function. Consider the set of lines through the centroid 
of G. It follows from continuity that at least one such line divides G into two 
convex subdomains of equal area over each of which the integral of w vanishes. 
We now divide each of these subdomains into two more convex subdomains of 
equal area over each of which the integral of w vanishes. 

Continuing this process, we arrive after a finite number of steps at a division 
of G into convex subdomains G, of arbitrarily small equal areas A,. Furthermore, 


(3.3) fudG= 


on each G,. 
Let 0, be the radius of the largest circle contained in G,. Then clearly 


(3.4) ALG, 
Hence, if A, is sufficiently small, the width o, of G is less than a preassigned e: 


(3.5) 0G: =. 


This means that G, is contained between two parallel lines at distance «. We 
introduce a rectangular coordinate system with the x-axis along one of these 
lines and the x,-axis tangent to one end of G,. Let L, be the length of the pro- 
jection of G, on the x-axis. Clearly, L,<D. Let p(y) be the length of the inter- 
section of G, with the line x,=y. Then #(y)<¢e. Because of the convexity of 
G,, p(y) is convex in y. 

Let M be a bound for the absolute values of w and its first and second deri- 
vatives. Then by the mean value theorem 


Ly 


(3.6) Sanbn is ) [w(0, ¥) Tay] <2M?A,¢, 
Gy ; 0 

(3.7) y) [w(0, »)] Pay] some, e, 

and 

(3.8) 


-SP0 )4(0,9)dy| SMA, e. 


Applying the inequality (2.10) of the lemma, we find, using (3.3) and L,<D, 


that 
y rf ou \2 
J leraduPace | aaa dG 
Coes a: 
20? D* [ wWdG —2M?(14+n®D*[14+ he) A é 


Gy 
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We sum these inequalities over all the subdomains G,. The sum of the A, 
is the area of G. Since « is arbitrarily small, we obtain the inequality 


(3.10) S| grad u|2dG = 22D fwd. 
G G 


Since # is any function with bounded second derivatives satisfying (3.2), we have, 
by definition, 
(3.11) We at LO 


4. The n-dimensional case 


Let G be a convex n-dimensional domain with boundary C (n=3). We 
again define w. as the infimum of the Rayleigh quotient (3.1) amomg functions 
having bounded second derivatives in G and satisfying the conditions (3.2)!. 

Let « be such a function. We consider the set of m —1-planes of the form 
aX, -4+6x,—=c passing through the centroid of G. By continuity we find that 
at least one of these planes divides G into two subdomains of equal n-volumes 
over each of which the integral of w vanishes. We divide these subdomains in 
the same way and continue the process until G is divided into subdomains G, 
of arbitrarily small n-volume V,. If 0, is the radius of the largest inscribed 
n-sphere, we have 


(4.1) V,1= K, of 


where K,, is the volume of the unit m-sphere. Hence, by a sufficiently large number 
of subdivisions, we can make g, less than a preassigned e. This means that each 
G, is contained between two parallel 7 — 1-planes at distance e«. In a particular 
G, we introduce new rectangular coordinates with the x-axis normal to these 
planes. We proceed to subdivide G, by means of planes of the formax,_,+0x,=c 
into subdomains on each of which the integral of « vanishes. We make these 
dividing planes pass through the centroids of the projections on x,—0 of the 
domains being divided. After a finite number of such divisions we obtain sub- 
domains G, whose projections on x,=0 have arbitrarily small »—1 volumes 
Vin-4_ It follows as before that if V"~" is sufficiently small, the projection 
of G, on x,=—0 lies between two parallel 1 — 2-planes at distance at most «. We 
keep the x,-direction fixed and choose the x-direction of a new rectangular 
coordinate system in G, perpendicular to these planes. 

If n> 3 we divide each G; further by means of planes ax,,_,+0%x,—=c passing 
through the centroids of the projections on x, = x,=0 of G; and of the succeeding 
domains. 

In this way we eventually obtain a subdivision of G into a finite number of 
convex subdomains G,’ over each of which the integral of uw vanishes. Further- 
more, each G,’ is contained in a parallelepiped of the form 


OS 7S, = IP Aten ed 


OS4,= 2, 


i 


(4.2) 


with respect to suitable rectangular coordinates. 


11Tf the boundary C is smooth so that the problem (1.1) possesses eigenvalues, 
[ty is the second eigenvalue of (1.1). 
Arch, rational Mech. Anal., Vol. 5 20 
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Let p(y) be the m—1 volume of the intersection of G,’ with x,=¥y. Then 
p(y) is convex because of the convexity of G,’, and p(y) Se"~* by (4.2). 


The inequality 
(4.3) fy = 7? D 


is now derived exactly as in § 3. 


This research was supported in part by the United States Air Force through the 
Air Force Office of Scientific Research of the Air Research and Development Command 
under Contract No. 49(638)-228. 
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Error Estimates for the Numerical Solution 
of Elliptic Differential Equations 


JOACHIM NITSCHE & JOHANNES C. C. NITSCHE 


Communicated by C. MULLER 


The numerical treatment of boundary-value problems for elliptic differential 
equations is often based on the difference method. The solution U produced 
by it, which is hoped to be an approximation to the exact solution w, is defined 
at the finitely many points of a grid and is obtained as the solution of an algebraic 
problem approximating the original boundary-value problem. The accuracy of 
this procedure will depend, of course, upon the fineness of the mesh. The esti- 
mation of the deviation U—vw in terms of this fineness is an important task. 
On the other hand, the quality of the error estimates obtainable finds its re- 
strictions in the regularity properties of the solution or, what is the same, in 
the behavior of the given data. 


In what follows the boundary-value problem 
Q (1s) = a(x, 9) Mey + 2B(%, 9) ey + 6(%, 9) My =F (x9) for O<%,y <1 


Vf 
”) u=0O for (x, y) on the boundary 
is considered. By replacing the derivatives u,,, u,, and u,, by difference 
quotients V,,, U,,, V., U;,, and V,, U,, of the grid function U={U, ,} in a square 
mesh of width 4=1/(N +1), asystem of linear equations L; ,(U)=£,,, is obtained. 
It is well known that an inequality |U—w|<CM,)h? can be derived, as was 
first shown, in the case b(x, y)=0, by S. GERSCHGORIN [3]. Here M, denotes 
a bound for the fourth derivatives of the solution uw. In order, however, that 
the solution u be four times differentiable and that such a bound exist at all, 
one has to assume that the coefficients a, b,c, and the right-hand side / possess 
continuous second derivatives at least. Unfortunately, these assumptions. are 
violated in many applications. Even in cases where a,b,c are well behaved 
the function /, often characterizing an exterior force, will not, in general, be 
continuous. Then an error estimate involving M, is useless. It is the purpose 
of this paper to derive error estimates under reduced assumptions. The coefficients 
a, b, c will be assumed to be continuous. From the continuity of the right-hand 
side f it can only be deduced that the second derivatives of u are square-integrable. 
But, as a matter of fact, the same statement is true even for square-integrable /. 
Therefore, / will be supposed to be an L?-function. Then for the values F;,, in 
the difference equations certain means of / must be taken. In cases where / 
corresponds to a point force, to a piecewise constant force, etc., these means 
can easily be computed. 

20% 
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The result, which in this introduction is given only for the case of a con- 
tinuous function /(x, y), is the following: 


(II) [Up — min] SC Lee (f) + IAI (ee (2 (6) + 2, (c) + #o*+0)]- 


Here the moduli of continuity ¢, are defined as «, (a) = Max |a(%, ¥1) —4(%2, Vo)| 
for 0%, Was Vip Vo teane aS %o|, [V1 — Va <p, etc. ||f|| denotes the L?-norm 
of f, and C is a numerical constant depending only on the constants of ellipticity 
of the differential operator &(w). For @ any number greater than / can be taken. 
With 9=/h' the last two terms in (II) contract to 2h*, and the bracket tends 
to zero as h—0. 

In the lemmas on which the error estimate (II) is based use is made of the 
theory of elliptic differential equations, especially of results of L. Bers [1], 
L. Bers & L. NIRENBERG [2], L. NIRENBERG [5], S. SOBOLEV [6]. The function 
spaces H,, introduced by P. Lax [4], prove to be a helpful tool. As a link con- 
necting the functions « and U, repeated mean values of u over squares with 
side 9=h (for instance g=h>) are introduced. Mean values over squares with 
side h have also been employed by H. WEINBERGER [7] for his study of K. O. 
FRIEDRICH’S symmetric hyperbolic systems. It seems, however, as indicated 
also by inequality (II), that the choice ga=h does not give the desired result. 


§ 1 
The boundary-value problem under consideration is 
(1) L(u) = a(x, Y) Uy, + 20 (X,Y) Mey + C(%, ¥) Myy =F(x,¥) for (x, y)ES 
“=O for (x, y) on the boundary dS. 


Here S denotes the square 0<%, y<1. The coefficients a, b,c and the right- 
hand side / are subjected to the following conditions: 


i) a, b, c are continuous in S= Saas: 

ii) There exist two constants g and Q (0<qSXQ) such that 
(2) qo + 9°) Sal + 2beyn+ cms QO(H + 7?) 

for all real numbers &, 7 and all (x, y) es: (uniform ellipticity). 

ill) f is square-integrable over S. 
Under these assumptions, as is proved in the theory of elliptic differential equa- 
tions, there exists a unique solution u(x, y) of (1). The second derivatives of 
are square-integrable functions, and the differential equation (1,) is satisfied 
almost everywhere. It is easy to derive estimates for the square norms of the 


second derivatives of wu in terms of the square norm of the right-hand side p 
In order to do this we observe first that the integral 


ff (Dee Pas v5) dxdy 


vanishes for every function v(x, y) which is of class C? in S and of class C2 in S 
with boundary values zero. This follows directly by partial integration from 
the identity v,, v,y— Ds = (V, Vyy)y — (Vy Vyy)y- Since uw and its first and second 
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derivatives can be approximated in the square mean by functions v with these 
properties and their derivatives, the same equality holds true for w: 


(3) JF (tas Myy — Way) dx dy = 0. 


Now aii eu the differential equation (1,) by w,,/c, we obtain 


(4) (aut, + 2bu Uy x Uzy + CU, y| + (u ca 4 Uy y —u,)=Lu,,. 


But 
Gat, 


(Ug + Usy) SOU t+ 2b ty, My + cour, S O(w,+u2,) 


on account of the uniform ellipticity. So, integrating equation (4), we find that 


(5) Off het dvdy ss ff ital dady. 


Introducing the ee! norms ||f||= (J ff Pdaxd y) etc. and utilizing the in- 


equality |« B| < — Sie apse ~~ B, valid for real numbers «, 8 and e>0, we can write 


(6) lee + LleeyllS SE Ueall® + 525 IAP. 
In a similar manner the other ne 
(7) lltesl2 + [eI +52 llAlr 


is obtained. If we set e=q?/Q and add the two inequalities (6) and (7) (omitting 
a term 2||w,,||?), we obtain, as announced, the estimate 


(8) Joell? + 2ll yl? + ly ll? 25> IAP. 


§ 2 
The system of difference equations is obtained in the following way: We 


cover S with a square mesh of width h= so obtaining a grid domain S’ 


with the boundary 0S”: 
SHEE SH Chk ih) AWS aR S N}, 


0S*:{P, = (th, kh);i=0,N +1 or k=0,N + 4}. 


1 
Net 


At all grid points P, of S” we define the difference quotients of a grid function 


= {U, ,}: 
Ver Uix = 55 (Uresn— 2U, 4+ Ua) 
(9) Vz Us,» = oa (igs, 041 — Usa, nga + Una, n—a — Gitar) 


V5) U; 5, =F (C4 20 a U; 4-1). 


296 JoacHim NitscHE & JOHANNES Cy GoNTTSCHE: 


The boundary-value problem (1) is replaced by the algebraic problem 
L;,(U") = 4; Vex U; pot 2 OV zy C5, + 6; 4 Vy yO, ,=F. for Peo, 


oo) U;,=0 for eS 


Here a;,, 0;,,,¢;,, stand for the values of the coefficients a,b,c in the grid 
points. The values of the function f(x, y) need not, of course, exist at the grid 
points P,,. So, if fis not continuous (in which case we can take F,,=/(ih,kh)), 
we shall choose for the F,,, certain mean values which will be discussed later. 

Our next problem is to derive a difference analogue of the a priori estimate (8). 
To this purpose we first observe that there exists a counterpart of equation (3) 
for any grid function V" vanishing on 0S”, namely 

N+1 7 


(14) (Vey Vi Y=S0, ) Vyy Vie) - 


4, k=1 1, 2! 


Vz denotes the unsymmetric mixed difference quotient 
(12) Vay Vin =z Vink — View ge Vy —1,k-17_ Vegas 
The difference quotient V,, V;,, can be expressed as a sum 


(13) Varnes =1(Vzy Vin + Vey V; +1; 1 ++1, pees Vi, eer 


SCHWARZ’S inequality gives 


~ 


(14) Vey Ka)* SPU, Via +, Vara (ey Maca a ee 
Combining (11) and (14), we obtain 


45) Df Cos Via) Woy Van) — Was Va 0- 
Now we proceed in the same way as before, obtaining 

; ay ia Vig Vasa) 20g View On wl Mey eel Pata eg Oy ae 
oe + {Pex U,,2) (Fy Ua) — Way T jek Oe 
It is useful to introduce the following norm for a grid function i 
a7) Ile = (ie 2 va) 
An argument similar to the one following equation (4) yields 
(18) Pee OU + 2 Poy Ol + Pro OES 2S FIR. 

§ 3 


Any real-valued L-function v(x, y) defined in S can be characterized by its 
Fourier coefficients 


(19) = Fa (x, y) sin(m x) sin(any)dxdy. 
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According to the growth behavior of the Fourier coefficients we introduce classes 
of functions: 


Definition. A function v is said to belong to class H,,, if the series 


(20) (1+ m? + n?)% (vo)? = |[o| [5 


i 


converges. 


Clearly, every L?-function belongs to class Hy, and its square norm coincides 
with the newly defined O-norm, so that we shall omit the index 0. 

If we denote by w”” the Fourier coefficients of the solution w of the boundary- 
value problem (1), then the inequality (8) can be written in the form 


foe) 
ely ee m4 (m? + n2)2 2 (umnye< i 


n=1 


or, using the fact that 1+m?+?< 3 a3 n*), 


[I/ll.- 


(21 ules 232 


In the following lemmas we shall derive a number of inequalities. For the 
sake of clarity we shall not always give explicit values for the numerical constants 
involved but rather denote them by C,, Cg, etc. 

oauue LN Lee paegEs T to class H, (x>1), then v is continuous in S and, 
for eS3, 


(22) Max | (x, y)| = C1() [lll 
(x,y) ES 
Cz (a) fipe is ai lle for a +2 
(23) é(v) S 


CxO (log + Nile Claw 10a ceo. 


Here «,(v) denotes the maximum of |v(%, 91) —2(%2, ye)| for |x— Xa\; 
lv. — ¥e| Se and (%1,%1), (%2, ¥2) ES. 2 


Proof. We first observe that the Fourier series >) v”” sin (a m x) sin(z 1 y) 
m,n=1 


converges absolutely and uniformly in S. This follows from 


a Jor S( 2d mens) Ct me + nny ome) 
m,n=1 m,n=1 m,n= 
r vay 
S20 f flo =~ 
0 
The same inequality yields (22) with C,( “=|. Further, again applying 


SCHWARz’S inequality, we see that 
|v(%4, V1) — U(%e, Yo) |? S o loll. 


=) (1+ m2 + n?)~* (sin (7 m x) sin (7 y,) — sin (7 m XQ) sin (0 yo)? 
m,n=1 
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Now (sin (av m x4) sin (a 1 y,) — Sin (2 7 XQ) sin (7 0 a) ye 
< 22 (1+ m? + n?) [Max (|x — *2l, 11 — Yel) )?- 
lence, for K 2 $, 
0687? [Max (|41— %ol, 41 — yal) 2 Mt + m+ n2)i—# + erie m> + n*)~% 


The first sum can be estimated as follows: 


R 
ay v dy 
Se er ca es (474 
m=+n?< R? Ps 
sory i+ Ry for «<2 
1+ R? 4(2 . y ot) 
=F [ “ds # log (1+ R®) for «=2 
1 
a ws ay ale veo. 
ie 
For the second sum we obtain 
2 PN ed ee vay 
ie ae Slt vies (1: 7?)™ 
42 = 42 A \o— 
od (1 + (R y2)?)? = a—1 (a) ' 


Assume that |x,—x,|<@ and |y,—¥,.| Se. In the case when «>2, let R->ov. 


Then as 


a eee 
_= ee 
In the case when «2, put Res , and make use of the inequalities 1+ R?<2 KR? 
and log (1+ R?)<4logR. Then ° 
@ cai ves for. G2 


5 
o ey n= er) Or << 2% 


The assertions of Lemma 1 follow from these inequalities. 


Lemma 2. Let v belong to class H, (x>1). Denote by vN (x, y) the finite sum 


(24) Oy > v™” sin(m x) sin(zn y), 


m,n=1 
and put W,,=0(ih, kh) —vN (ih, kh). Then 
(25) I| Mlle S Ca (@) 2 Ilo fla. 


Proof. The grid function W,, can be expressed in terms of its “discrete” 
Fourier coefficients 


N 
(26) wn" — 4h? >" W,, sin (th mi) sin(ahnk) 
eS 
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in the form 


N 
Wr= d W"" sin (ah mi) sin(xhnk), 
m,n=1 
and 
N 
(27) 7 Wile=a 2 (en)! 


Using the orthogonality relations 


i si for w=m(2(N+1)), 
(28) > sin (a h m1) sin (0 h wt) 


= a —5, for w=—m(2(N +1)), 
0) for w= +m(2(N +1)), 


which are valid for 1S mSN and all w=1, we see that 


(29) We — >)” 
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Here v"” are the Fourier coefficients of v(x, y), and the sum has to be extended 
over all w= +m(2(N+1)), y= +n(2(N+1)) with the exclusion of the pair 


“=m, v=n. Then, again applying SCHWARz’s inequality, we see that 


N N 
SD ees 0 > Ot wt S(t ptt ote (oe? 
m,n=1 m,n=1 pw=+tm w=tm 
(14,») (mn, n) (11,») +(m, n) 


1l<m,nsN 


< ol, Max 5 Whitey 


v=itn 


(14, ¥) = (m, 2) 


Now, for any m and in 1Sm,n<N we have 


tarts D+ +1) e+ 9), 


pein b+q2l 


ve4%n 
loo) (oe) 1 
ia 2 pa ip ean (8 Pod) ef 


(4, ») + (m, 2) 
se dy 
<i) [+f Ji) a a ‘et, 


‘ meta fhe ea 
af 9 
=a +, — 
Hence, with (27), 
(30) Ws (5 +n) P lle. 


arend: 


, 
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Lemma 3. Let v belong to class H,, for an integer «22. Then the second der- 
VAlIVES Vyyx, Vyy, Vyy 4S well as the second difference quotients’ V,, v, Vy 8 Wage 
belong to class H,_,, and 


(31) lleeellema>--+7llFog Plla—a = 7 [IP |la.- 


Proof. The statement follows immediately for v,, and V,,v (and v,, and 
V,,v). Namely, the Fourier coefficients belonging to these functions are 


— mv" and Tele sea 1d 
ahm \? 
a 
The mixed derivative and difference quotient can be represented, in the square 
norm, by series Re 
>) mnv™" cos(m x) cos (mn), 
m,n=1 


me s sin(ahm) sin(ahn) 


mn 
SiON hn ee cos (7% m x) cos(an y). 


m,n=1 
Since the norm ||v,, ||, 2 does not change its value if we use the Fourier coefficients 
with respect to the system {cos ( mm x) COS(mn y)}, it follows that 


Il? ylla—2 = 26) m? a2 (A+ m2 +e (oma)? << lly 2; 
m,n=1 
the estimate for eee follows similarly. This proves the lemma. 
For later purposes we subject our functions to a smoothing process. We 
define a mean 
v 


a+ 
(32) ogo) (x9) =se ff v(m) dé dy. 


Q 


In order that this mean be defined for all points of S, we continue the function 
v across the sides of S periodically as an odd function. The relation between 
the Fourier coefficients of v and those of «, (v) is 


TLC ONS IE92:0 
sin 2 sin 2 
(33) Wen O) ate yan 
2 2, 


oye Se ae or 
(34) [le (2) lla S|] 2 lla 


Lemma 4. Let v belong to class H, (c=0). Then W,(v) belongs to class H,4,, 
and 


(35) le lla S375 Jol. 


* They are defined, similarly to (9), by Kexu(*, y) = 5 [u(~+h, y)—2u (x, y)+ 


v(*—h, y)] ete. 


Error Estimates for Elliptic Equations 301 


Proof. On account of the inequality m?n?=4(1+m?+ mn?) we have 


(1+ m+ m2) (00,(v)"™")* <5. (o™)? 
and hence (35). q.e. d. 
Lemma 5. Let v belong to class H, (1<a<3). Then, for oS<#, 


; 2(a—1) 
a Max. |o(#,y) — eo) (1, 9)] Colo) oo 


Proof. We have 


v(x, y) —@,(v) (%, 9) = > Smn? v™” sin (mm x) sin (zn y) 


m,n=1 
where 
sin eee Se @ 
i D 2 
Smn il + mm n 0? 
Hence 
2 
eae Min {2, we tm +n}. 

Then 


|v (x, ¥) — @,(r) (x, »)| S 4] le ul  ghalt +m? + ?)~* 


m,n=1 


If R= 4, the sum can be estimated as follows: 


Dy Bon (it 0? + 02) Ty ai 
m,n=1 +nv2sR* m®+n*?>R? 
< oe ‘5 (1+ m? + n?)?-% + 4 yy (1+ m? + n?)—% 
576 m+n?s R? m+n? >R? 
Now 


R+y2 
(4 1+ m? + n?)?—- “Ss ac AS ois ae ee Gee ee 
+s 


m+ n?<R?* 0 


fepsie zie Min (2, a) Ct aia 


Further 


[o-e) C—L 


m+n? >R? R-y2 1R 


Consequently 


co 2 T0n 2 (3—Mi (x, 2) Tt 
> Sant 2 +) a mare Vase DS 


Bee 


m,n=1 
2 
Putting R=o ax), we obtain 
4(a—1) ae On 
2 2) —a < 9 Max (2, 
D shalt tm 98) ae rg eee eg 
mM, n= 


which proves (36). 
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§ 4 
The solution u of the boundary-value problem (1) belongs to class Hy. Its 
second derivatives and difference quotients belong to class Hy, and in general 
the corresponding Fourier series will not converge everywhere, in particular not 
necessarily at the grid points. This difficulty can be overcome by applying the 
smoothing process of the previous section. We introduce the function 


(37) ii (x, y) = 09 (w) (%, ¥) = Oy (o, (wt) (%, ¥)- 
ti (x, y) belongs to class H, according to Lemma 4, and 
(38) Il@|laS Coo |IFI| 
o- 
6 ng? V2 


Lemma 6. Let D® be any second partial derivative and V? the corresponding 
difference quotient. Then 


(39) ||D?@ — V? ||, 5 Ch? o*||f\I. 


The grid norm in (39) is taken for the grid function D? “(th, kh) —V? u(th, Rh). 
N is defined as in Lemma 2. 
Proof. According to Lemma 2 we have (a= 2) 


|| D2 a — D2 GN ||, < Cy (2) || D2 all, 
72a — V2a" ||, <Cy(2) WIIDP@ 


ls, 


and hence, by Lemma 3, 


(40) |D2a— Dar ||,,  ||V?%# — VW ||, SC, h? 9-4 |[7|| 
(CRON) ec C,). Further 
N 
Dea — par P<2 d 2, (am 
m,n=1 
Here 
4 sin (72) 
702 m2 \1 hE IB Peas hou 
sin(ahm) sin(ah ? 
Lnn= mm n (4 ( “)) if)" =40. 
4 sine (74) 
72 n?\ 4 mh nt if D* = d,, 


For all these g,,,, the inequality 


2 
3 MUL m? + nd)? 


Emn = 


holds true. Therefore 
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or, using (38), 
(41) || D2 a = V2ar||,S me h2 


=< UC, 2-4 
eo lall: 
Combining the inequalities (40) and (41), we obtain 
||D? a — 724i], <||D°a — Deal, + ||DPa — 72a", + ||2a — V2", 
rm AG, a 
Se ||7I|(2C, +4) =Celo [l/l 


that is, our assertion (39). 

Let a(x, y) be a continuous function in S. Denote by e,(a) its module of 
continuity: ¢,(@)=Max |a(%,, ¥;) —a(%2, ¥2)|, where the maximum is to be 
taken over all pairs of points (%, y,), (%2, y2)€S with |x,—x2|, |v, —y2| St. 
It is easy to extend a(x, y) into the whole (x, y)-plane in such a way that the 
module of continuity is not increased. In the following lemma it will be assumed 
that a(x, y) is so extended. 


Lemma 7. Let a(x, y) be continuous, and let v belong to class Hy. Then 
(42) || cog (v) — (av) ||, S 42, (2) |||]. 


Proof. We have 


{ame (v) — w3(av)} (x, ¥) 


@ 

2 

f 48 dn! {a(x,9) AC ASUS a ecu ie. 
2 xu(x été, ytntn’). 


Applying ScHWARz’s inequality to the inner integral, we obtain 


[fa (8 (v) — o8(»)} (x, | Sea)|ta ff en) asan 


7 +e yt+e | 
e— 059 —0 


Hence, using the fact that e<([o/h]+1)h, for e=A, 


2008 (2) — 08 (a 0) |b SM (e_ (a)? 45 4([ 2] + 4) iinet as 


S 16 (¢,(a))? I ‘ 
This proves Lemma 7. 
After these preparations we can turn to our estimation of the error. We have 


(43) Lx (@) = [2 (@) Jeineny + RY, 

with 

(44) R= {4,4 (Vex — On2) + 20;a Vey — Oxy) +. 6,2 Voy — yy) } BA, RH). 
But, on account of the definition (37) of 4%, 


(45) [2 (%) |cin,en) ee [os (A) lin,an ate RP 
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where 


(46) RP, = f(a? (thy x) — 03 (A Myx) + 2b (05 (Hey) — 9 (6 Mey) + 


+ (609 (yy) — 9 (C Uy y)) Sink h)* 
According to Lemmas 6 and 7 and the inequality (8) (without omission), we get 


(47) [RP UlpS Max (|a], [0], |el) {Il Pex— 22) Ale + 21] Pry — 2) Mle + I(Vyy— yy) Hla 


[FR lpS 4 eq (4, B, 0) {| 4ec|] + 2lleyll + Iles ll} 


(48) f 
< 40% (a, 0,0) [Al 


Here ¢,(a, b, c)= Max (¢, (4), & (2), &(c))- 

We compare the solution {U,,} of the difference equations (10) with the 
values of the function % at the grid points, because their difference satisfies the 
linear equations 
(49) L,,(@—U) =a (f) (th, kh) —Fin + Ret Rea. 

Consequently, by virtue of inequality (18), we obtain the main inequality 
I|P.2 (% — U)|le + 21, xy ( (hie U)lle + [Vay wha U) z 
< Cll?) — FB + Calif (€e(2, 8, 6)? + Cu 0 [L/P 
60? 6Q4 96 04 C3 
C= 8, Cp = , Cy = an 

The final exploitation of this main inequality (50) depends upon the properties 
of the function f(x, y). We distinguish two cases: If /(«, y) is continuous in S, 
then the numerical calculation is carried out with the values f;,=/(vh, kh) for 


the right-hand side of the difference equations. Using Lemma 7 with a=/, v=1, 
we then find 


(50) 


Here 


I] eo(/) —F 


5 44,(f). 
Thus we obtain 


MN ee ee Mie + I[Vyy(# — O) Ike 


54 
OH SIF? (Co (a, b, c))® + Cy ht o~ 8) + Cra (& (7)? 


(Cia 16.C,): 
In case f is not continuous, we must keep the expression }||@3(/) —F |], in 


place of ¢,(f). Of course, if for the numerical calculation we use the mean values 
we(f) (¢h, kh) for the right-hand side F,, then this term does not appear. 


§ 5 
The main inequality (50) of the last paragraph is now used to derive the 
estimate for the difference u(ih, kh) —U,,. First we prove 


Lemma 8”. Let W,, be a grid function with vanishing boundary values. Then 


(52) Max |Wiel S2V1+2||Vee+Vyy) We. 


1Si,kSN 


* This is a finite difference analogy to an inequality of S. SoBoLEv. 
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Proof. As usual, we express W,, in terms of its discrete Fourier coefficients 


N 
W..= >»; w"" sin(a@hmt)sin(ahnk). 
Mm, n==1 
Then we find 
ad : ahm ; ahn\\2 
I| Vex + Vyy) W le = pape oe +sint(#>*)) (Wns: 


On the other hand 


pmale( 3 perc St at) 5: owt mye ores 


m,n=1 


1 ~ 1 
SMVex+ Vy) a Da aman 


since 


M 
z 
+/— 
SS 
II\ 
|= 
N}aQ 
— 
— 
Le 
z 
a) 
a 


m,n=1 i 
This proves Lemma 8. 
The difference u (th, kh) — U;,, can be written as 
[uw — We () |cin, kn) a [w, (u) — ws (1) lin, an at [a (th, kh) i U; 4). 
Applying Lemma 5 for «=2 and using (34), we obtain 


Max _|u — a, (u)| <C;(2 ) o||x|le, 
(«, y) ES 
Max |, (u) — w9(#)| SCs (2) @ ||, (#)|l2S Cs (2) @||%llo- 


(4,9) €S 
Therefore, Lemma 8 in conjunction with the main inequality (50) gives 


Max |u(ih, kh) —U,,| 
(53) 151,kSN 
<Cygl|09(f) — F lle + || fl] {Cra & (4, 8, ¢) + Cis? o™* + Cre of. 
Here 


Cye=_ 20a) Cy, Cu= P2042) Co, 


Fie EN Cl 2D 
C= = V20+ 2) Cy, C15 = ae C5(2). 


In particular, if we set =A’, 


Max |u(ih, kh) —U,,| 


(5 ) 1<1,kSN 


<= C3 l|one (/) 14 €y3 (4, 9, ¢) ial + (Cis + Cis) hi Wale 


2(f) if the function f is 


As before, the term ||@%3(f) — 
continuous. 


This paper was prepared under Contract Nonr-710(16) between the University 
of Minnesota and the Office of Naval Research. 
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Feblerabschatzung 
fir die numerische Berechnung von Integralen, 


die Losungen elliptischer Differentialgleichungen enthalten 


JOHANNES C. C. NITSCHE & JOACHIM NITSCHE 


Vorgelegt von C. MULLER 


Vorgegeben sei das Randwertproblem einer elliptischen Differentialgleichung 


© (u) = a(X, Y) yy + 20(%, Y) Uy + C(%,¥) Uyy=F(x%,¥) in S:{0<2x,y<14} 


4) 
(1) u—O aufdem Rande @S. 


Zur numerischen Integration wird (1) durch ein System von Differenzenglei- 


chungen fiir eine in-den Gitterpunkten eines quadratischen Netzes mit der 

7 ; definierte Gitterfunktion U’={U,,} ersetzt. In den 
= , 

Anwendungen treten haufig Integrale der Form 


(2) ee y) u(x, y)dxdy,..., bios Y) Uyx(%, y)adxdy,... 


Maschenweite ) = 


auf. Bei der numerischen Berechnung wird man diese Integrale durch Summen 


N N 
(3) Wi rides fi ge Ee) ASE NEIA GY Fe 
ersetzen. Es ist das Ziel der vorliegenden Arbeit, Abschatzungen fiir die dabei 
resultierenden Fehler zu geben. Eine Abschatzung fiir die Differenz U; , —u (th, kh) 
ist bereits in einer vorangehenden Arbeit! hergeleitet worden. Dabei war Wert 
darauf gelegt, die Voraussetzungen fiir die Koeffizienten a, b,c und die rechte 
Seite / — welche wir hier iibernehmen — méglichst gering zu halten: a, b und c 
werden in S=S+4@S stetig und f in S quadratisch integrierbar angenommen. 
Die Differentialgleichung soll gleichmaBig elliptisch sein. 

Im Einklang damit wird von der Funktion #(%, y) vorausgesetzt, daB sie in 
S quadratisch integrierbar ist. Im allgemeinen braucht dann # in den Gitter- 
punkten nicht definiert zu sein. Fiir die in den Summen (3) auftretenden Werte 


1 Error Estimates for the Numerical Solution of Elliptic Differential Equations, 
vorangehend in dieser Zeitschrift; im folgenden mit ¥¢ abgekiirzt. 
Arch. rational Mech. Anal., Vol. 5 m1 
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P, , werden vielmehr gewisse Mittelwerte eingesetzt, namlich 


(41) Ba=or)in= LLL AOR S ie 


2 

Der Parameter 9 kann dabei geeignet gewahlt werden. Falls p(x, y) stetig ist, 
wird man natiirlich 
(4.2) Eon=b (th, kh) 
nehmen. Wir werden die Fehlerabschatzung erst fiir den Fall (4.1) durchfithren 
(Ungleichung (30)). In § 5 wird dann die auf (4.2) beruhende modifizierte Fehler- 
abschatzung angegeben (Ungleichung (32)). Wir beschaftigen uns, was offen- 
sichtlich geniigt, nur mit den Integralen tiber zweite Ableitungen von w. 

Nehmen wir fiir einen Augenblick an, daB und f stetig sind, so 1aBt sich 
das Resultat folgendermafen ausdriicken: 


SS p(x, y) Dudxdy — 1S plik, kM V?U;, » 
iS i,k=1 
SC [DIL eg (f) FIA & (B) + &9 (B) eo (Ff) HINDIII AII (€(4, 8.) +A? 0-4). 


D?u ist eine zweite Ableitung von uw und V?U der entsprechende zweite Diffe- 
renzenquotient der Gitterfunktion U”. || || bedeutet die L?-Norm, und die «, 
sind die Stetigkeitsmoduln. Indem man o=h” (y<#) setzt, strebt die rechte 
Seite von (5) mit / nach Null. Falls alle auftretenden Funktionen einer Lipschitz- 
Bedingung geniigen und man y=2 annimmt, so la4Bt sich der Fehler durch die 
Schranke Ch* abschatzen. 

In den nachfolgenden Ausfiihrungen werden die Bezeichnungen der Arbeit 
Jt beibehalten und die dort gewonnenen Resultate benutzt. 


(5) 


§1 
Um die spateren Entwicklungen nicht unterbrechen zu miissen, soll der 
folgende Hilfssatz vorausgeschickt werden. 


Hilfssatz. Die Funktion v(x, y) set in S quadratisch integrierbar, und es 
werde V(x, y)=we(v) (x,y) gesetzt (9h). Dann gilt 


(6) IV eS 4 Mell. 
Beweis. Es ist 
Q @ 2 
v |, 2 2 
V Ce fp) PE d d , , , , 
Ive % > elf Edy [[ de dif oh +E+e, kh tytn) 
; irs ite 
2 2 


Die Zahl of =h([ 2) =} 1) ist gréBer als 0, und folglich 
- : Q 
JJ ai dy vGh+é+E, kh+y +m Soff didn Pih+é&, kh+n'). 
@ = 


2 
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Damit wird 
Iles mer > S [fever PUh+E, Rht+n) <4ie (2) + 4) [olf 


2@ 


Es ist aber +1 >» Woraus Ungleichung (6) folgt. 


& 
h 


§ 2 
Zur Herleitung der Abschatzung (30) bedienen wir uns wie frither der aus der 
Lésung «(x, y) des Randwertproblemes durch zweimalige snwendnng des Glat- 


tungsprozesses w, gewonnenen Funktion %(x, y). Der Ausdruck h? > PV? U; 


1,k> 
welcher sich anter Einfithrung des skalaren Gitterproduktes rat 
(7) Us W)= 3 ViaMs 


abkiirzend in der Form (P, V?U), schreiben la8t, wird in der folgenden Weise 
umgeformt: 


(8) (P, V2 U), = (P,V?(U — a)), + (P, (V2 — D*) ui), + (P, D*i),. 
Zufolge Lemma 6 von % gilt 
(V2 — D®) a||,S Ce 


so daB 
(9) |(P, V? —D?) &),| S || Pll, ||(V? — D?) H||,S 4Cg h? o-* |||] |[FII- 


Bei der letzten Ungleichung ist der in §1 bewiesene Hilfssatz benutzt. 


Fir das erste Glied auf der rechten Seite von (8) erhalten wir mit Hilfe der 
Hauptungleichung (t 50) 


[(P,V7(U —%)),| ||P lle {VCs Ileoe (Z) —F ||, +V Cro Il fl] & (a,b,c) + Cy, W? o- ale 
bzw. mit y, = 4 Max ices Vor, GHP 
(10) |(P,V?(U —%),.| Sy [lA ll {llee(?) — Fle + [l/l] (o(@, 8,6) +? o4)}. 

$3 


Der dritte Ausdruck rechts in (8) laBt sich in ein Integral verwandeln. Wir 
fiihren den Kern 


(11) P(x, 93) =4 Dd Sn Sin (7% x) sin (7m y) sin (x m &) sin (7 0 7) 
m,n=1 

ein, wobei 

oe woes 

sin 
(12) Smn — 8m bn» =( av 0 

2 
Dann ist 
(13) &(E,n) = SS y(n x, ) u(x, y) dxdy 
S 


2c 
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und 
(14) (P, Dea =I OZR P, a Yo (ih, Bhs x, y)} (Du) dx dy. 
Wir werden nun zeigen, daB die Funktion 
N 
(15) TG, I) TE EASY) 


fiir N->oo im Quadratmittel gegen die Funktion £(%, y) konvergiert. Auch hier 
gehen wir wieder schrittweise vor. Indem wir 


N 
(16) ph’ (x, 93 & n) =4 D Sm, Sin (xm x) sin (70 y) sin (zm €) sin (a nn) 
m,n=1 
setzen, definieren wir 
N 
(17) Fi(4,y) = WD Ba Ye (oh, bh; x, 9). 


Dann gilt 
T(x, y) — T(x, y) ‘ 
= 4)" on, Sin (mm x) sin (any) ) i" UR 
k= 


m,n 


*e Sin (7 hm 2) sin (wh nk}. 
1 


Hier, wie auch im folgenden, bedeutet 


N co N loo) oo 
i Sa By SNS 
= M 
m,n m=1 nm=N+1 n=1 m=N+1 mn=N+1 


Weiterhin ist 
[a= al? a (T(x, y) — T(x, y))2dxdy 


P N 
= 42) ban [l? vie 


Unter Beriicksichtigung der Periodizitatseigenschaften des Ausdruckes in der 
geschweiften Klammer und der Abschatzung 


+, in ( (ah mi) sin (hn k)} 
al 


Sm < Min (1, iam) 


m2 (a m2 


ergibt sich 


1 
aie tall 
yg S 2 
= ‘ Gin Sant q(N+1) ey (—1)49* P,sin(athmt) sin (hn k)| au 
n= 1k= 


2 
(— 1)?" P., sin (hm) sin (xhna)h = 
1 
2 2 = 2 
§m+-p(N+1) 8a+q(N+1) > (—1)?*+9* P sin (hmi)sin (xhna)} 
1 i,k=1 


N 
1/74 <a pa") | Saet Gee ec 
sz (#2 Ph) |e a + (SES Ay 
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und folglich 
TAZ a 
(19) 7 — 5) <2 eo |p). 


Nunmehr beachten wir, daB P., mit den Werten der Funktion w;(p) in den 
Gitterpunkten (7h, kA) tibereinstimmt: 


(20) E> San he simithm 1) sin(annk). 
n=1 


Dabei steht g,,,, fiir die in (12) erklarten Ausdriicke, und #”” sind die Fourier- 
Koeffizienten der Funktion #(x, y). Wir schneiden die Reihe (20) ab und defi- 
nieren die weitere Gitterfunktion 


N 
(21) PX, = Dd Sinn B"” sin (7h mi) sin (x hn k) 
R=1 


m, 


und entsprechend eine Funktion 


(22) TA %s y) = HD BN yf (Ch, hs, 9). 


Pit 
Die Differenz der Funktionen J, und 73 ist 
N 
T, (x, y) — Tp (x we ie Linn Pr” sin (ahmt1) sin (hn k) x 
N 
a 4D Sy Sin (er hei) sin echo 8) sin (x x) sin (av y) 


N , 
rat Dy Gay Sin (oe pt Xx) sin (zv y) » SPER payt pe ° 
il 


Ly mM=rpE 
n= 
(m,n) + (H, ») 


Die Kongruenzen sind modulo 2(N+1) zu nehmen, vgl. (% 28), (9¢ 29). 


N 2 
ee Cay) 2. SE eee 
Bv=l ME=—re 
n=+y9 
(m,n) + (u, ») 


Anwenden der Schwarzschen Ungleichung und Beriicksichtigung von g,, <1 
ergeben 


1syp,vSN 


= 
i | 
S 


|Z—nIPst | > s “(omy 
(m, n) = (u, ») 


Das Maximum 148t sich ahnlich wie in (18) durch (= h2 0 z: abschitzen, so 
daB wir 


Lg Snap cos ns 
(23) |. — Bis 42 0 [IAI 


erhalten. 
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Driicken wir in dem Ausdruck (22) fiir Tj die PN und Yo durch ihre trigono- 
metrischen Darstellungen aus, so erhalten wir 


N 
=> 2 ,p""sin(am x) sin(an y). 
m,n=1 


Daraus ist ersichtlich, daB 7; (x, y) nichts anderes als die abgeschnittene Fourier- 
Entwicklung der Funktion 
=> 2pm" sin(am x) sin(an y) 
™m, n=1 


ist. Somit ergibt sich 
|| — a (0)|P= a 2" Smn (b™")?, 


und da im Summationsgebiet der Summe a der rechten Seite 


Tp oe 
ey =( i 0 ‘)' 
ist, weiterhin 


(24) |Z, — 4 (6) || < 29 # oI]. 


Als letzte Abschatzung bendétigen wir noch die fiir die Abweichung der Funktion 
we(p) von p. Es ist 


io — of P= 2d (1 chal? 


Unter Verwendung von 


1 — gan < Min (1, © (om? + n°) 
ergibt sich 


|e — ws (p)|?<— a 0 4 (1—A) Ds (pm) + + ys (pmn)2 


m?+n?<o-24 m24n®>o-2A 
(0<A<1). Wir fithren die Abkiirzung 
(25) dalPl =e eee al 
m+n >R 


ein. Fiir jede L?-Funktion # konvergiert 6g(p) nach Null, wenn R nach Unend- 
lich strebt. Nunmehr wird 


(26) Iba II <2 


7 GA ||p||+6,-2(8)  (0<A<4). 


§ 4 
Wir fassen jetzt die vorangehenden Abschatzungen zusammen. Aus (8), (9), 
(10) folgt 


|(P, VU), — (P, D?a),| 


(27) 
es [loo (A) — F lle + IA ILIA 1 €0 (4, 8, ¢) + 72h? on} 
(Yo=%1+4C,). A 
(P, D®ii), — an (x,y) Deudxdy 
(28) 


= JI {(T—-T) + G-T) + (L—o8 (6) + (w§ () —p)} (D2 u) dxdy 


S) 
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finden wir mit Hilfe von (19,) (23), (24), (26) und den a priori Abschatzungen 
(Xt 6, 7) fiir || D?«|| weiterhin 
eae — LF ble, y) Dudxdy| 


<2 7 (2 mero 28 rele += elloll+ 4-s(0)} 


Dabei haben wir A = } gesetzt. 


Indem wir die Ungleichungen (27) und (29) kombinieren, erhalten wir schlieB- 
lich die gewiinschte Fehlerabschatzung 


bal (x, ¥) Oy bynes 


ee 
Sy llelllos) —Flle+ aoe 


+ ||P ll Ill fe, (@ 8 ©) + ys? 0-4 + 4.0} 


“ 2240 , 160 _ 
(ys=724 45 q? Pigeges Oe rl 


§ 5 

Die Fehlerabschatzung (30) ist unter der Annahme, daB die Funktion p(x, y) 
quadratisch integrierbar ist, abgeleitet worden. Dabei hatten wir fiir die Gitter- 
funktion {P,,} die Werte der zweimal geglatteten Funktion # in den Gitter- 
punkten genommen. Falls # eine stetige Funktion ist, dann wird man natiirlich 
P.,,= (th, kh) wahlen, so da die Formel (30) erganzt werden mu8. Wir wollen 
abschlieBend die Fehlerabschatzung fiir diesen Fall angeben. 

Die bei der Herleitung der Ungleichung (29) verwendete Abschatzung (26) 
der Norm von w3(p) — p kann bei stetigem # durch den Stetigkeitsmodul (vgl. 2, 
Lemma 1) ausgedriickt werden. Zundachst ist klar, daB der GlattungsprozeB den 
Stetigkeitsmodul nicht verschlechtert: 


&, (co (D)) ss e(B) 9 =A, 2,20) 
Weiterhin ergibt sich 


_Folglich haben wir 
el) — Pl] Max |oc(?) — pl 24, (2). 


(ey €5 
Lemma 7 von % liefert (im Falle v = 1) 
05 (6) — Plle S 4&0 (6 
Somit wird unter Benutzung des Hilfssatzes von § 1 
|(P, V2 U), — (p, V? U),| = |(we() — 2, V? U),| 
(1) S 4e,(6) ||V? Ul S <P (0) | 


2 


I 


ge cal + Ge eel) — FDU 
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Kombination der Ungleichungen (30) und (31) gibt 


N 
ILS p(x, 9) Deudxdy —18Y p(th, kh) V2U,, 
S 


i,k=1 


(2) < (yllAl| + <2 ¢(6)) 7-8) lle + +P Ill (6) + 
+ [PII AIl ty &0 (4, 8, ©) + 75h? a4}. 
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1. Introduction 


The literature concerning unsteady motion of a body in a fluid with free 
streamlines can be divided roughly into three categories. 

The first category deals with exact solutions. In that case it is necessary 
to restrict the field to unsteady motions of a special type. VON KARMAN [7] 
was the first to take that direction. Later publications appeared by GILBARG [2] 
and Woops [3]. 

In the second category we shall class those articles in which the unsteady 
motion is considered as a perturbation of a steady flow. In this way Woops [4] 
investigated the unsteady motion of a flat plate according to the Helmholtz- 
Kirchhoff model. In [5] Wu gave a formulation of a boundary-value problem 
for the Roshko model (also named transition model). Closely connected with 
this is the research on the stability of free streamline flows. ABLow & HayEs [6] 
published a report on the subject. Their results were extended by Fox & Morcan 
in [7]. 

The third category contains papers in which the whole flow is linearized. 
TuLtn [8] was the first to apply the method of linearization; he considered a 
steady flow. PARKIN [9] solved the unsteady problem of an infinite cavity 
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behind an inclined flat plate and gave a formulation of the finite cavity case. 
In [10] and [5] Wu developed a theory for the unsteady finite cavity case. 
Timman [11] and the author [72] treated the case of a partially cavitated inclined 
flat plate. The unsteady motion was considered small with respect to the steady 
flow (second linearization). In order to get a closed cavity it was necessary to 
admit a discontinuity in the unsteady vertical perturbation velocity in the 
wake of the body. Another possibility 1 is to drop the condition that the cavity 
be closed in the unsteady case and to require that the normal component of the 
unsteady perturbation velocity be continuous in the wake. The last possibility 
will be investigated in this paper by comparing the linearized cavity flow with 
the corresponding reentrant-jet flow. For a steady cavity flow it has been 
proved that the linearized theory yields a first order approximation to reentrant- 
jet flow theory [13]. From the following it will appear that also in the unsteady 
case a linearized theory can be formulated in such a way as to become a first- 
order approximation to the reentrant-jet theory. We shall prove this connection, 
just as in [13], for the simplest case, that of a flat plate moving perpendicularly 
to its plane. 

The paper is divided into two parts. In the first part we consider a non- 
linear steady reentrant-jet flow. The unsteady motion is treated by first-order 
perturbation theory. Since in this theory an arbitrary motion can be considered 
as a superposition of harmonic oscillations, only a harmonic time-dependence 
is introduced. 

In the second part we treat the linearized problem. The solution is obtained 
by two different methods. The first uses the complex velocity and is directly 
connected with the method of the first part. The second works with the complex 
acceleration-potential and offers the best prospect for more complicated problems. 

Finally, it should be remarked that in this introduction we have omitted 
papers that deal with the initial stage of impulsive unsteady motion. Several 
Russian authors have worked in that direction, e.g., GUREVICH [/4]. Their 
investigations are directed towards the determination of induced mass of the 
cavity (see BiRKHOFF & ZARANTONELLO [15]). In this paper, however, we are 
interested in the unsteady motion as a whole. 


Part I. Reentrant-jet flow with unsteady perturbations 


2. Formulation of the problem 


A flat plate is placed in a flow of an incompressible non-viscous fluid with 
uniform velocity perpendicular to the plate. Behind the plate a cavity with 
reentrant jet has formed. The jet is supposed to disappear on a second sheet 
of the physical plane considered as a Riemann surface. This artifice is used 
to represent within the framework of potential theory the dissipation of energy 
that is observed at the rear end and in the interior of the cavity (see BIRKHOFF 
& ZARANTONELLO [15]). Here we prefer the reentrant-jet model to the Ria- 
bouchinsky model, since the former shows more resemblance to experimental 
observations (see EISENBERG [16}). When the plate performs an unsteady 


* This possibility emerged from a discussion with Dr. P. EIsENBERG and Dr. M. P. 
TULIN. 
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motion, the shape and size of the cavity will change. If we assume that the 
deviation from the steady position is small in a certain sense, this unsteady 
motion can be considered as a perturbation of the original steady flow (lineari- 
zation of the unsteady effects). 


First we treat the stationary ie @) 
flow, and then we investigate the Ser ea act 
effect of the unsteady motion. Vee Bes) 
Y, 2 
3. Steady flow - Beer, | 
The steady symmetric flow 2 pena Seacrest i 
about a perpendicular flat plate Fig. 1. Reentrant-jet flow about a flat plate 


with reentrant jet is shown in 
Figure 1. The Cartesian coordinate system has been chosen in such a way 
that the middle of the plate lies at the origin, and the velocity of the fluid at 
infinity has the direction of the positive x-axis. We denote half the breadth 
of the plate by b. The fluid is assumed to be incompressible and non-viscous, 
and the motion, steady and irrotational. The relation between the pressure and 
velocity is then given by the Bernoullian law: 
Pot toU=p+seP=b,+20%, (3.1) 
where g = magnitude of velocity, 
p = pressure, 
U = velocity at infinity, 
poo = pressure at infinity, 
9, = velocity on the cavity, 
, = constant pressure on the cavity, 
o = constant density of the fluid. 
We choose the units of length and time in such a way that the constant fluid 
velocity g, on the cavity is equal to 1. The susceptibility of the fluid for the 
occurrence of cavitation is characterized by the cavitation number o: 


Gen he (3.2) 
+0 U2 
From (3.1) and (3.2) it follows that 
U= (1-40)7?. (3.3) 


The x and y-components wu and v of the velocity vector g satisfy the conditions 
of incompressibility and irrotationality : 


ae - =), 
Xx iM 

4 
oy au aa 
Ox oy 


These equations show that w= —iv is a holomorphic function of z=x+7y 
in the interior of the fluid. On account of (3.4) we can define a velocity potential 
gy and a stream function y such that 


ss Ge 
Ox oye” 3 5) 
oy oy ; 


a = 
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Then ®=+7y is a holomorphic function of z=x-+7 y, and 
w= d® ; (3.6) 


For the moment we consider only the upper half of the physical plane, 7.e., the 
upper half of the region occupied by the fluid. As a result of the presence of the 
reentrant jet this region is two-sheeted. Since the region is simply connected, 
Ir @ we can map it conformally onto the 
4 upper half of a ¢-plane ((=o+7T) ?. 

| The images of three points on the 
| boundary can be fixed arbitrarily. We 
| 


A B C 1 g Choose the image of z=0o as t=on, the 
(-a,) 0 7c,0) (do) +image of B as t=O and the image of 
Fig. 2. Conformal image on the upper half of the ¢-plane A as t=—a, where a is a fixed positive 


number (see Figure 2). 
In paragraphs a) and b) below we determine ® and w as a function of ¢. 
The conformal mapping giving z as a function of ¢ is then given by 


dz=wid®=w1 as, 


or 
t 
i a®@ 
— =1 oe 
z= | w rate (3.7) 
—a 
: : 4 : LD aD. 
a) Determination of ® as a function of t. We consider rap ta ile holo- 


morphic function of ¢ in the upper half-plane. Since the real axis is the image 
of streamlines (where y =constant), oe is rea] along the real axis. By means 
of the reflection principle of SCHWARZ we can continue ae analytically into 
the lower half-plane. a can have singularities on the real axis. These cannot 


be essential singularities (for a proof see BIRKHOFF & ZARANTONELLO [15]). 
Singularities can occur only at t=co and t=c. From a consideration of the 


local behavior it follows that 4 remains finite at too and behaves like —© 
6. 
at ¢=c. As a result of the branching of the streamlines at the point D, ae 
will have a simple zero at t=d. Thus 
d®D > t—d 
dh Tey 2) 
where K is a real constant; hence 
@ = K|t — (d — c) log (t — c)] + constant. (3.9) 


b) Determination of w as a function of t. We introduce Q=log w=log g—79%, 
where # is the angle between the velocity vector and the positive x-axis. Q is 


» Of course this o is different from the o defined by (3.2). This will cause no 
confusion. 
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a holomorphic function of ¢ in the upper half-plane, where it satisfies the follow- 
ing mixed boundary conditions: 


—o<a<-—a, 9=0 Oia lint == Or 

+ 020, G=5n -or ImQ-=—4n, 
OO <6. G4 or RKeQ= 0, (3.10) 
€<6 <4, G=—7 or ImQ=z, 
ax fax ce) v= Or.) Lin — 0, 


For the solution of this boundary-value problem we continue 2 analytically 
into the lower half-plane by means of the reflection principle of SCHWARz: 


2() = 20. (3.11) 


The boundary-value problem is now equivalent to the following Hilbert 
problem (see MUSKHELISHVILI [/7]), viz, to determine Q(t), holomorpnic in the 
upper and lower half-planes, satisfying the following linear relations on the real 
axis: 

—coco<—a,’ 2 —Q-=0, 


—a<o<0, 2Q*—-Q=— x21, 
Vie aC, Grae — Ov (3.12) 
BO Kitt QO* — O-=2n1, 
d<0-< 00, Bi 2 = 0. 


Here 2* and 2 denote the boundary values that are assumed by 2 in approach- 
ing the real axis from the upper and lower sides. Moreover, 2 has to be continuous 
at the points ¢=0 and t=c. At infinity 2 must remain finite. 

We consider first the homogeneous Hilbert problem, 7.e., the above-mentioned 
problem with the right-hand sides in (3.12) replaced by 0. A solution Q, of 
this homogeneous problem is 


2,= jit —°). (3.13) 


The cut imposed so as to make the square root single-valued is taken along 
the real axis from ¢=0 to t=c. The branch is chosen in such a way that 


Vit—e)~t as t>oo. 


The Hilbert problem (3.12) can now be formulated as follows: Determine Q(é) 
in such a way as to satisfy on the real axis the relations 


Ori 1% 
hes Bhai Sek te as | iS ot 
Onis Or rt 
ae O, == 
Hage ar We V2 a(Ge) 
Ons Qi 
ieee eel kad 0, (3.14) 
shonturelg eae Se Qni 
Qj On| 2 Veta 
Colak em (OM at Sp 
Acero aos) ha 
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By means of the formulas of PLEMELJ (see MUSKHELISHVILI [17]) we find that 


seg Ie _ =(¢—a) (0-2) tot f att ae eal ae 


After some calculation we arrive at 


(Vi(e+ alt )] 1 [Va (¢—c)—Vt(d—c) | IG 
eg OE Tare a = VA Se alee ea ae ieee ee 
and : ‘ 
[Ve(c+a)—Va(t—c) ]* [Va(t—c)—Vit(d—c) ] a7 
5 ‘(Viet ae c) 1? [Vd(t—c)+Ve(d—c)] O) 
The expression for z as a function of ¢ can now be written as 
atk Wicca cele [Va(¢—c) + Vt(d—c) ] t—d rah 3.18) 


[Vé(c+a)— c) 1? [Va(t—c)—Vi(d—c)] *—¢ 


Remark. (3.16) and (3.17) represent 2 and w also in the lower half of the 
physical plane when we map it conformally onto the lower half of the ¢-plane 
by means of (3.18). 

In the solution, which is the combination of (3.17) and (3.18), three unknown 
constants, viz, K, c and d still occur. These unknowns are determined by the 
following conditions: 

1. at t= co we must have w= U=(1+0)7, 

2. z has to be a single-valued function of ¢, 

3. the distance between A and B=b= half the breadth of the flat plate. 

We first formulate Conditions 2 and 3 mathematically. The condition that 
z should be a single-valued function of ¢ reads 


eat (3.19) 


Cc 


where C is a closed contour surrounding the plate and cavity, 7.e., surrounding 
the part of the real axis that lies between —a and d. Application of the theorem 
of residues to the exterior of C gives 


{res ene a 0. (3.20) 
Condition 3 can be written as 
0 
he ee (3.21) 


We now apply the three conditions to the solution obtained. Condition 4 gives us 


ue ei Va} Wa— Vd—c] 


4 
( aaa ~ (Vera+Val* (Vd+Va—c] : 


(3.22) 
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Condition 2 furnishes the equation 


14/- cld 
Seale eta laa fe C= (3.23) 


From Condition 3 we obtain 


° (Val (ca) ee (o— c) ]8 [Vd(o—c) Vo(d c)] o 


} (Vo( (c+a)— Va(c—c) |* : [Vd (o—c)—Vo(d—c) | Os 


“do=bi. (3.24) 


Elimination of a between (3.22) and (3.23) gives the following relation 


ye+Voeth-Veyesy. ox 


Moreover it can be shown that 


een ae ee 
d—c=28 (14 a0 \j2 ire], (3.27) 


The integral in (3.24) can be expressed in terms of elementary functions, but 
since the result is rather unwieldy, it is omitted. (3.25) gives the relation between 
c and o. From (3.27) follows the connection between d and o. Finally (3.24) 
determines K as a function of b and o. Thus the solution for the steady case 
is completely determined. 


between c anda: 


1+o= 


and 


4. Unsteady motion as a perturbation of steady flow 

Let the flat plate perform an oscillation of the form he’®’ normal to the 
main flow. By T we denote the time variable; / is the amplitude of the oscillation. 
The deviation is taken positive in the direction of the negative x-axis. Thus 
the normal velocity in the direction of the negative x-axis is equal to jwhe’®’. 

Remark. It is tacitly assumed that in all expressions containing exponential 
terms complex with respect to 7, we must take the real part in order to get the 
true time behavior. In the formulas this real part is denoted by “‘re’’. Since 
we linearize the unsteady part of the motion, this procedure will cause no trouble. 
Of course the imaginary units 7 and 7 are unconnected; we have 7?=—1 and 
j2= —1, but 77 is not reducible. The real part with respect to 7 is denoted by 
‘Re’, the imaginary part is denoted by “‘Im”. 

We now suppose that wh is small with respect to g, and that / is small with 
respect to 6. In this case the unsteady motion may be considered as a pertur- 
bation of the steady flow. In order to formulate these conditions more rigor- 


ously we introduce two non-dimensional parameters ¢, and é,: 


f= a (velocity parameter), (4.1) 


and 
aya & (amplitude parameter). (4.2) 
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The conditions mentioned above can be expressed by 
SR iid ae WS ag OS (4.3) 


To express the relative magnitudes of ¢, and ¢, we introduce the parameter 


== te (frequency parameter) 
9 Ve 


and require that 
Pao (4.4) 


The reason for this requirement will be discussed below under a). We shall 
neglect, as usual, second and higher powers in the disturbances in comparison 
_ with first order terms, and moreover terms proportional 
Ze to €) in comparison with those proportional to ¢,. This 
B last procedure is motivated by (4.4). 
Juhesol The boundary conditions will now be formulated. 


a) On the flat plate the boundary condition is 
(see Figure 3) 


“ , fi 
1 : ohe® 
Gus 2g Bares = 


5 (4.5) 


We put P=PYot Pi, P= Pot Pr, I=Ptrn, V=I+A 

and p=)+/,, where quantities referring to the steady 

7 flow are denoted by the suffix 0, and those referring 

Fig. 3. Unsteady velocityonthe tO the first order terms of the unsteady perturbation by 
flat plate the suffix 1. These first order terms include, for the time 

being, the terms which are linear with respect to ¢, or é9. 


In a first order approximation (4.5) may be replaced by 


c h pier 
j= — Oy , 4.6 
= Z— 8,4 2% (4.6) 
This approximation is not valid near the stagnation point, but nevertheless, 
as is customary in linearized theories, we assume it to be valid there. We wish 
to apply this condition to the unperturbed position of the flat plate. In a first 
order approximation it is true that on the flat plate 
a8 


De een 
Do 5 a h, (4.7) 
where the derivative “2 must be taken at the unperturbed position of the 


plate. Using the fact that log gj —7®% is a holomorphic function of z, we may 
reduce (4.7) as follows: 


o = He 0 log qo pee 0 log qo h 

2 oy 2 OPo 8 
4. 
Tt Fo p, — Tt Ody 1 O-—=6) o (4-8) 
2 OM 2 60 K o-—d, * 


In this reduction use has been made of (3.8). 
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Remark. cy and dy denote the steady values of c and d. This distinction is 
necessary, since as a result of the unsteady motion c and d will change with time. 


Boundary condition (4.6) may now be written as 


. jo T 
} _ johe! | Odo 1 Oo— Co 
on aaa. (4.9) 


In accordance with condition (4.4) the second term in (4.9) may be neglected 
in comparison to the first term at the right-hand side. Thus the final form of 
the boundary condition on the flat plate is 

0, = jah ed?™ . 


x (4.10) 


The condition must be applied at the unperturbed position of the plate, 7.e., 
on the interval —a<o<0O of the real axis of the ¢-plane; g,) as a function of 
o is given by 


z = (4.14) 


% [Va (cy 0) +V —o(cy+ a) }* (Vd, (cs—c) + V—o(d,—) | 


_ Wa(ce,—o)—V—a(e,+ a) ]* [Vay (650) —V—o (a—e) | 


This follows from (3.17). In the sequel we shall write Q(c) for Ge 
0 

We shall now discuss briefly the difficulties encountered in the case when 

€, and ¢é, are of the same order of magnitude. Then the second term on the right- 

hand side of (4.9) must be retained. Since in the neighborhood of o=0 the 


expansion of g,, as derived from (4.11), begins as 


“ty 4por Cota yoo 
ss ]/ a Co +2) 


the behavior of ae there is seen from the following term: 


rai ea | Cota ae 1 At 
=p glee yargenw te 4-13) 


Thus becomes infinite at the edge of the flat plate. This fact can be for- 
oO 


mulated physically by the statement that there the acceleration of the fluid is 
infinite. The occurrence of this singularity 1s known in the potential theory 
of steady cavities separating from sharp edges of obstacles. This singular be- 
havior, however, is not realistic for actual fluids because of the viscosity. Now 
solving the boundary-value problem by taking (4.9) instead of (4.10) as the 
boundary condition on the flat plate, we can see from a formula similar to (4.40) 
at the end of this section that the perturbation pressure #, is no longer con- 
tinuous at the edge of the plate. From physical reasons, however, we should 
expect continuity there. Therefore it is plausible to conjecture that in unsteady 
motion when ¢, and e, are of the same order of magnitude, the point where the 
cavity separates from the plate is no longer fixed at the sharp edge of the plate. 
Since we wish to restrict our attention to cavity flows with fixed separation 
points, we have introduced the assumption (4.4). 
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poems, (4.12) 
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b) On the boundary of the cavity we consider the Bernoullian law for unsteady 


flow: 


se + SP tel ac(n. (4.14) 


Splitting up all quantities into steady and unsteady parts, we may replace 
(4.14) to a first order of approximation by 


Se Stee + 4+ = C(7), (4.15) 


Since S ga + ae = constant Merten: the entire region occupied by the fluid, 


we may reduce (4.15) to 


se +g + B= DT). (4.16) 


If a =r,, the Bernoullian law reads 
i) 
o +an+2 =D D(T). (4.17) 


Remark. In partial differentiation with respect to T the space variables 
are kept constant. This is equivalent to keeping q, and wy constant. 


On the boundary of the cavity we have #,=0. Thus 
0 
o% + gy = D(T). (4.18) 


We wish to apply this equation on the boundary of the unperturbed cavity. 
Letting points correspond to each other by orthogonal projection, we find that 
(4.18) may be replaced to a first order of approximation by 


re +%=D(T), valid on the unperturbed cavity, where gg=1. (4.19) 
On the oe of the unperturbed cavity, however, 
Wo = const. 


Differentiating (4.19) with respect to @y, we get 


64 OD meee 
OT Od, | 0% 


ee (4.20) 
We prove first that 
4. 7%, —1%, is a holomorphic function of ¢, and 


OD, 
— =. 2? 
stan, (4.21) 
2. Putting 6B= @,4+ @,, we have 
d®,  d® UO We _ w 
d®, 4d, Pree dP, a Wo igs (4.22) 


To a first order of approximation 


we 2 __ pv +2, = 6% (44 01) =v, (1+2,), 


or 
——1=Q. (4.23) 
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Combination of (4.22) and (4.23) gives 
Q sere, (4.24) 
Always to the same order of approximation, 


Q = log w= log — 19 = log {go(1+ )}— 184 — vd, 


(4.25) 
= log gg ty —10) —18,=2,+%—-1%,. 
Since 2 and Q, are holomorphic functions of ¢, we conclude that 
2,=71-18, (4.26) 
is a holomorphic function of ¢. From (4.24) and (4.26) follows finally that 
ne = Re QO ==7;. q.e.d. (4.27) 


On the unperturbed cavity the boundary condition is now 


On, OV ae 
OT ra Lio ee (4.28) 
or 
7, =G(IT — @); (4.29) 


1.¢., 7, is.a function of T —q,. 
In the case of harmonic time-dependence we get 


r= gel), (4.30) 
where g is a still unknown constant, real with respect to 7, but possibly complex 
with respect to 7. 

In considering here this harmonic time-dependence, we anticipate, since justi- 
fication follows only by (5.15), expressing g as a linear function of wh. 

For Q=2,+, the following mixed boundary-value problem can now be 
formulated: Determine (¢), holomorphic in the upper half-plane, in such a 
way that on the real axis it assumes the following boundary values: 


—o— oe — a4, F=0 or ImQ=0, 

—a< 6 <0; O=41ntjwhQ(c)d?? or IMQ=4n—jwhQ(o) e?", 
OZ60 <6, log gq tn= grr or ReQ=g e%lF—m), (4.31) 
CRA ORS Gd: 0=—2 or Jim =a 
C2 0= oj} DEO or ImQ=0. 


Putting Q(t) =Q(0) according to the Schwarzian reflection principle, we can 
transform the mixed boundary-value problem into a Hilbert problem with the 
following linear relations along the real axis: 


=o o—— 4, 07 .0R-- 0, 


4.2 0e0,, OF = O> = 71 = 217 ah O(c) f°, 
O<o< Cc, Or + Q- = 2¢ lot), (4.32) 
Eee teal he OF OS 77, 


6 =acoy we hea >), 
22* 
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Remark. To a first order of approximation g can be taken equal to zero at 
Cy<a<c if cy<c. Putting c=—cy+c, and d=d)+d,, we can split this Hilbert 
problem (4.32) into the following Hilbert problems for 2, and Q,, where 


Q= 2, FQ; 
~w<cx—a, 23—-2,=0, 
—a<o<0,  2—-2,=—x1, 

O< 6 =e, 25+25=0, 
C=O 2G 2, — 2, = 271, 


d<¢ <te, V = 2, =0, 


and 
= © ag = Oo, te ee, 


—a<o<0O, Qt —Q; =— 277 wh Q (a) e'?7, 


OL OA kee ante eae 
Cp SOE co, Oi Oe, 


According to (3.16), the solution of (4.33) is 


[Va (t—c) —Vi(d— 


(Vd (t—c) + Vi(a— 


(4.33) 
(4.34) 
c) J 
Ta (4.35) 


We now determine the solution of (4.34). By means of the solution of the homo- 
geneous problem 2,= Vi(t—% t—c») the Hilbert problem (4.34) can be written as 


follows: 
ae [2,)* 1.2. 
PRM on ArT 
‘te. ee us | Q, es are é 
a<a<o, | o*| —|g¢] = — 24/0206) 
OF apr = KOR 2g ef? (T—) 
Ox 0 Ce. ‘| ‘| = 
; | Qh 2, QR 7 
“sot le abe 


When 0<o0<Cy we have 


Po = K [a — (dy — Cy) log (cy — )]. 
Q (a) is given by 


(4.36) 


(4.37) 


Q(c) — Wale oe o(cota)]” | [do(eo—2) +V= alo] (4 gg 


‘[Va(cy ey pas (Cg + +a) |* [Vay (Co o) V a (dy 


Co) ] 


Using PLEMELJ’s formulas bs MUSKHELISHVILI [/7]), we get finally 


Q(o)do 
Co—9) (o—2) 


Q, = ae eT Vi(t =e) if ae 
a ( 


S. 1 OL 1/57 é 
—£¢ yee) f 
0 


—j w K[o—(dy)—cy) log (co—o)] 


Vo(¢)— a) (o—t?) 


do. 


(4.39) 
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Substituting cy+c, for c and d,)+d, for d, expanding and combining the ex- 
pression obtained with (4.39), we find 


Oc a1/ t >| a +// dy ze Vi(@—cy) di a3 
Co V t—eyo|2 V cote dy—Co d, t—d, dy—Co 


0 
j@h joTigis .\ Q(c) do 
4 JO4 oT VE —¢ (i 4.40 
e Ke sh V—o(e9—) (o—2) Say) 
x? —jw K[o—(dy—Cq) log (co—s)] 


eae Wi(een ye 
ae Vet Co) i do. 


: Vo(cy—s) (o—2) 


5. Additional conditions 


In the expression (4.40) for 2, three unknown quantities still occur, viz, g, 
c, and d,; g is a constant, possibly complex with respect to 7, since we have 
already anticipated the time behavior of G(T —qp); the two other unknowns, 
c, and d,, are functions of 7, complex with respect to 7, which will appear to 
be also of exponential type. All three quantities are real with respect to 7. For 
their determination we have the following three additional conditions. 


1. The perturbation velocity must vanish at infinity, 7.e., 
Q7=0 as |t] oo. (5.1) 


This condition gives us, as a result of the symmetry of the problem, only one 
condition real with respect to 7. 

2. The perturbation cannot have a source or vortex at infinity. The first 
part of this condition is equivalent to the requirement that the pressure should 
remain finite at infinity, which generalizes the closure condition to the unsteady 
case, for a source at infinity with time-dependent strength gives rise to an 
infinite pressure there, according to (4.17). The second part is trivially satisfied 
as a result of the symmetry of the problem. Therefore this second additional 
condition gives us, in the same way as condition (5.1), only one condition real 
with respect to 7. As can be seen from (4.24), it is formulated mathematically by 


resQ2,=0 as |t|—> oo. (5.2) 


3. Finally we have to ensure that the pressure on the cavity is constant not 
only in space but also in time. This last condition is not expressed by the boundary 
condition, since the boundary condition (4.28) has been obtained by partial 
differentiation with respect to gy. We shall now formulate Condition 3 mathe- 
matically. 

The Bernoullian law for unsteady flow has been written in the form (4.17): 


OP 2 i 
pie tee ae ee (5.3) 
We set g,=0 when |t|/—>co. We are allowed to do so, since , is determined 


only to within a function of time and q, is finite at infinity as a result of the 
absence of a source there. Since at infinity 7, =0 and ~,=0, we have D(T)=0. 
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From (4.24) we infer that 


n= Re [2 to,—Re [4h team Re [2,4 st dt. (5.4) 
On the cavity the following identities are valid: 
I=1, r= geerF—e), i, = 0. 
Thus on the cavity (5.3) can be written as 
901g folT—) (5.5) 


oT 
or, by use of (5.4) and (4.37), as 


oO 
K Re [ O25 ys fy g i®{T—Klo~(4—) 04 (@—9)}}, where OSa<cy. (5.6) 
MMe PSG, 

—oco 
This is the mathematical form of the third additional condition. It is easy to 
prove that this condition is independent of the choice of o in the above-mentioned 
interval. 

Since the three additional conditions are linear in the unknowns, and since 
the known members of the equations depend harmonically upon the time, we 
may put 

q=6,d°?, d=d,de?. 
For G(T —q@p) we have anticipated this time behavior as already remarked. 
€, and d, are real with respect to 7 but in general complex with respect to 7. 


If we let C=O e/°™ (5.6) goes over into 


o 
Rek [Q,7= di= J ge joke or ohoelerel = je gece Stew 
—oo 
The path of integration may be deformed according to CAUCHY’s theorem. 
We shall now use the three additional conditions in order to determine the 


three unknown constants, ¢,, d, and g. Applying (5.1) to the expression (4.40) 
for Q,, we get 


fi E aoe | Sa a, |) do 
Cal al eeecea Ue VM dee. dy V do—hq 
0 


iw h a (oj K[o— (dy—c5) log (co—s)} pe 
jo eS 5 0 — Co) log (cg—o 
oe eee ae CO). 
It sh e( ) V—a( 6 (cy—0) ath Va (cy—o) : 
Applying (5.2) to (4.40) and using (5.8), we arrive at 
’ b j= et gt 
“e) or 2: d aa see (5.9) 


se fo odo 4 & tf j@K[o—(dy—cy) log (¢)—2)] odo 
=9) Vo (cy—0) 
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From (5.8) and (5.9) we show immediately that 


A 


41 a do dy 
B Vets oe £0 1) 
(5.10) 
ae foo Coach A hi ph oe fe fenle a atentes oN] e212? 
6 (¢,—0) ms Vo(cg—o) ’ 
AVF 
= 
; aes i (5.11) 
=i" f Ee Ef erie Rtete—ados tao) Co 8) Eo 
nae ¥=e(¢)—0) ay Vo(cy—a) 


By means of (5.10) and (5.11) the expression for 2, can be reduced as follows: 


(Om we t(t (do—9) (Co— 9) oh 
Leis) af ow Vo (¢.—2) (a—#) (tcp) (tay) ae 


(5.12) 


jo K(o9~(d—c)108(6-9)] (gg) (ga) 
ew o do. 


ie 
pay San fee 
) ( “df Vo (cy—a) (o—t) (ty) (¢—dp) 


The first integral in (5.12) can be expressed in terms of elementary functions, 
but since this result is not needed, it is not given here. We are now able to apply 
Condition 3 as expressed by formula (5.7). We first take o=0 in (5.7) and use 
formula (5.12) to substitute for OF The left-hand member of (5.7) is now re- 
placed by a repeated integral. Since the integrand is a holomorphic function 
of ¢ except at the end point ¢=0, the integration with respect to ¢ can be per- 
formed along a path which meets the interval [—a, cy] of the real axis only 
at the point ¢=0. We now move the end point of the path of integration from 
t=0 to a neighboring point t=e, with Ime>0. Since both integrals in the 
repeated integral are now absolutely convergent, we can change the order of 
integration. Evaluating the inner integral and afterwards letting ¢—0, which 
is legitimate on account of the uniform convergence of the outer integral, we get 


ieee 2 0 dy—o o— fo g Lco=e | 3G" Fe eed 
ie Clo beac eo Sierra ary 
= ie 2 Peet gine cae bl 
: : (5.43) 
ey 20 dy—o 
—j @ K[o— (dcp) log (¢o—3)] ‘ im 
+ fe |are cos| a 1) ae do 
0 
=i IE gi@K (doer) log co 
@M 
Use has been made of the following identities: 
Relig Vi(t—cy) dt en, =O 4 ie Veo — =a! if <0, 


ona (ty)? (6-4) Cy by OBO Wesel) —=o] 
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= 2. — [are cos( 1)||/ z “ tf 0Se05cy 15-1924) 
Co—o & Co—5 Co—O 


Remark. This expression is finite for o=cg. The last integral in (5.13) can 
be reduced as follows: 


and 


£0 


TE ai eo (¢o—9)] ATG cos ( 20 — 1)| dy—o do 
Co Cy— 9 
0 Co 
= a larc cos (2° = 1)| d e710 la— (do~e0) log (co—4)} (5.14) 
oO 0 


—joK[a—(d.—e,) log (¢o—a)] 
= aay, ei okla 0 — C0) log ey me hope 
Lg o = a(¢)—9) 


The integral occurring in the last member of (5.14) and the integral of the same 
type in (5.13) can be expressed in confluent hypergeometric functions. 


Using (5.14), we can write (5.13) finally as 


—jon| [ov dy—o 2-1 tn [0104 log Weo—o o— Veal, ao| ! 


Ve-—c+V—o 
Ss (5.15) 


tafe ef eKlo— (dy—Cy) log (¢o—«) ] Ay gu. gu 2(d, 
Vo(cy—o) 


0) Jig|do=0. 


Formula (5.15) gives us g; ¢, and dy are then found from (5.10) and (5.11). The 
solution is now completely determined. 


6. First order approximation for small breadth 


«) In this section “‘small” is understood to mean small with respect to the 
length of the cavity. We start by determining the first-order approximation 
for the steady flow. For that purpose we fix the value of K at 1 and consider 
the quantity a, which has so far been an arbitrary constant, as an unknown 
to be determined from the condition (3.24) of Section 3. Combination of this 
formula (3.24) with (3.25) and (3.26) or (3.27) gives us the relation between a 
on the one hand and 0 and o on the other hand. For the computation it is con- 
venient to consider a and ¢ as given and to determine 0} and o as functions of 
a and c. We suppose a to be small. It will appear then that for fixed c-0 also 
b is small; z.e., as a tends to zero, b also tends to zero. The cavitation number 
o as a function of a and ¢ is given by (3.25). We expand this expression in 


powers of \2 so that 
14 om 1+ 4 [Spo (6.1) 


Hence we have to a first order of approximation 


o=4 \/2 : (6.2) 
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We now expand the integral in (3.24) in powers of \/2 , using (3.26) and (3.27). 


For that purpose we introduce first a new variable of integration and take the 
real part of the integral. 


[Va (c—o)+V—a(c+a)]* [Vd(e—o) + V—a(d—o)] o dag 
[Va(c 0) V a(c a) |* (Va (c—o)—V de) | O-6 (6.3) 
1 2 
my ee a(d—c) 1—y|” 4ac(c+a) (1—y?) 
=| [1+ SSS ee) ee ay, 


where the new variable of integration y is given by 


= Va(c—o)—V—o(c+a) 
Va(c o)4 V o(c+a) : 


(6.3 a) 


According to (3.26) and (3.27) we have to a first order of approximation 


G6 6 A a 
in C 


and 
de =) (6.4) 
Hence to a second order of approximation the expression in (6.3) becomes equal to 
1 
alee aaa ae 
4a f yt dy =F (m+ 4). (6.5) 
0 


From this formula it is seen immediately that to the same order of approxi- 
mation 


b=—a(n+4). (6.6) 


Expanding the expression (3.17) for w into ascending powers of \z. we get 


‘a || fie i 
Sac ca <a ee 
Substitution of (6.2) gives us 
o ||/t—< t 
a a / ria Va 


These expansions are uniformly valid in the interior of the cut f-plane. It may 
easily be seen that to a first order of approximation 


ee, (6.8) 


di— a1, (6.9) 


Hence to that approximation w is given by 


= | Ve} (6.10) 
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The relation between 0, c and a is 


o=4 eae (6.11) 
In the neighborhood of z=0 the behavior of w is given by 
w= VE Lt eraiee (6.12) 


The quantity c may be considered as a first-order approximation to the cavity 
length. 

f) First-order approximation for the unsteady motion. We determine the first- 
order approximation to the unsteady perturbation field for small breadth of the 


plate by expanding (oh in a series of ascending powers of /z. From (5.12), 


by again introducing y as a new variable of integration, we find that to a first 
order of approximation 


1 
ies jok ea as, 
; Nie ARIE Se 


: eV ol: eat on cae (6.13) 


= lah ee é mi 
f ma Vea i =a. 


- e122 ae — 
ae = 


Relation (5.15) connecting wh and g is equivalent, to the same order of approxi- 
mation, to 


. — a 2-9) — = 
—2oh(n- 2 |) Cte¢g fe = de =U. (6.14) 
R 


Formula (5.10) becomes 


= a £o > Sty 
os | a qaoh a 4 —jwo 
ig fa 2st eta Sa e fm EE tor ony 
0 


The integrals in (6.14) and (6.15) can be expressed in terms of cylindrical func- 
tions by means of the following identities: 


joo 4g 


[igen Flo) 649 


and 


Co ee 


fom ata ny i lo) ih (o“)). (6.17) 


0 
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Hence (6.14) may be written as 


—2joht2 1/4 + ge F[7,( 2) (1 — L) +) h(o 2}] =0 (6.18) 


@ Co 


and (6.15) as 


“A ea 5 Co 
4/4 =2j0h2 2 egies 3), 


Lo Co T% Co 


From (6.18) and (6.19) we show immediately that 


g=2jah2t2 / aM J° 2 hyo “2 (4 i bil (02) (6.20) 


Tt +4) Cy @ Co 


and 


eet at2 lax Co fi . Cy\|—1 
OC, =2wh = Jo(w 2) | Jo (o ) (4 jh(o 2) = WG(624) 
In the derivation use has been made of (6.6) in order to eliminate a. For wd, 
we find the same expression as for w¢,. The relation between w¢, and g reads 


3 x 2b me Facts Co 
jo |epyy Wee” Blo 4) 
or (6.22) 


We shall now discuss the conditions to be imposed on the parameters w, h 
and 0, in order that the approximations introduced so far be valid. 

Remembering the discussion at the beginning of Section 4 and maintaining 
the notation introduced there, we must first require that 


ae eee (6.23) 
oe (6.24) 
a= et. (6.25) 


As a result of the approximation for small breadth 20 of the flat plate, we must 
further require that 


a= e<t. (6.26) 


However, to be sure that the unsteady effects are not of the same order of 
magnitude as the second-order steady effects, we must impose the additional 
condition that 

Mie (6.27) 


£9 Ic 


or, in terms of the smallness parameters ¢, and ég, that 


4X &. (6.28) 
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The five conditions (6.23) to (6.27) are not independent. An independent set 
of conditions is given by 


Rss (6.23) 
Came (6.25) 
ee Sars, (6.28) 


Using (6.6), we may replace (6.28) by 


= 2b 6.2 
a= | enn <* Oa 


From these conditions it can be shown easily that 
fi (7 wai! \ o 1 6.30 
max 4, Sah ta a ae (6.30) 


heh ZG; (6.31) 


In the case that the approximation for small breadth is not carried through, 
the conditions (6.30) and (6.31) have to be replaced by 


and 


1 @ 1 
Sea ke (6.32) 
and 
h<b. (6.33) 


Returning to formula (6.21) and taking into account that according to (6.30) 
w>>1, we conclude that the velocity of the moving rear end of the cavity is 
nearly in phase with the velocity of the flat plate. The deviation of the rear 
end of the cavity from its steady position is of the same order of magnitude as 
the deviation of the flat plate. 


Part II. Linearized theory of unsteady cavity flow 
7. Summary 


In the linearized theory to be discussed in the following sections, the vertical 
flat plate is replaced by a singularity. It will appear that the velocity field of 
this singularity in the steady case is identical with the field given by (6.10), 
1.e., by the first-order term of (3.17). The strength of the singularity can then 
be expressed in terms of half the breadth 6 of the flat plate. 

In the unsteady case the strength of the singularity is assumed to vary with 
time. Considering the variations in this strength to be small compared to the 
steady part of it (second linearization), we may proceed in a manner analogous 
to that for the steady case. It will appear that in the case of harmonic time- 
dependence with circular frequency m the unsteady part of the flow field can 
be identified with the first order term of On as derived in Section 68, if w is 
assumed to be large and the amplitude of the unsteady part of the strength 
of the singularity is put equal to 
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A treatment of the same problem by means of the acceleration-potential 
requires the application of additional conditions, different from the ones used 
for the method of the complex velocity. It is shown in Section 11 that both 
methods are equivalent, provided the right conditions are applied. 


Finally we give an expression for the drag in the case when our linearizing 
assumptions are valid. 


8. Linearized theory of steady flow 


We suppose that the flow pattern as a result of the presence of the flat plate 
can be considered as a small perturbation of the original uniform flow field. 
This will generally be the case if the breadth 20 of the plate is small with respect 
to the cavity length. The plate can 
then be replaced by a singularity. le 
The rear end of the cavity will show | (Z) 

a similar singularity. Within the 
frame work of linearized theory the | 
occurrence of these singularities | 

| 


can be inferred from the presence — _v=0_ u=0 _W=0_  & 
of a stagnation point in the neigh- (6,0) (1,0) 
borhood. We shall not go into the Fig. 4. Linearized boundary-value problem 


details but refer to [1S]. 


We choose the Cartesian coordinate system in such a way that the cavity 
extends between the points (0,0) and (/,0). The velocity at infinity has the 
direction of the positive x-axis. The units of length and time are so chosen 
that the constant velocity g, along the cavity is equal to 1. We put g=(1+-4, v). 
Since 1— U<11, we have u<1 and v<1. We shall neglect second and higher 
powers of w and v and their products (linearization). w=«—z71v is a holomorphic 
function of z=x+7y. 


The linearized Bernoullian law reads 


PP. a Ch (8.1) 
Hence we show that 
Ug = — FO. (8.2) 
On the cavity, 
u=0. (8.3) 


Since the maximum width of the cavity is of the same order of magnitude as 
the perturbation velocities, we may apply condition (8.3) on the x-axis between 
(0, 0) and (J, 0) instead of on the cavity. The holomorphic function w now must 
satisfy the following mixed boundary-value problem (see Figure 4): Determine 
w(z), holomorphic in the upper half of the z-plane, such that it assumes the 
following boundary values on the real axis: 


See a 1 () ea) sm Ore. Line (), 
O24 0,” ~or Rew=0; (8.4) 


=a co st), Of lmw—0. 
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By means of the Schwarzian reflection principle this boundary-value problem 
for w(z) can be transformed into a Hilbert problem with the following linear 


relations along the real axis: 
=O 40, 7 = Y= 0, 
O<x<l, wr+u=0, (8.5) 
LG AO: — ae Oe 


The general solution with integrable singularities at z=0 and z=/ is 


preted (8.6) 


se z—l 


The cut for both square roots is made from 0 to /; the branches are so chosen 
that 


and 


V5 ~1 as [z| > oo: (8.7) 
We assume that the strength of the singularity at z=0, determined by the 


flat plate, is given. There we have 


wr Ali. (8.8) 


=e 


Hence the strength of the singularity =A yi. This gives us one relation between 
A and /. For the determination of the unknowns 4, B and 7 we still need two 
relations. These are furnished by the following additional conditions: 


1 w=— 570 at z=, according#o(6.2), (8.9) 
2. rés.W=0 at z=00. (8.10) 


This condition is equivalent to the requirement that the cavity be closed (see 
[15]). From (8.9) we deduce 


A+B=-—fto, (8.11) 
and from (8.10) 

A—B=0. (8.12) 
Hence 

A=B=—}o. (8.43) 


Thus for w we have the formula 


ae a9 


This expression agrees with (6.10) if we identify c and 7. From (6.9) it follows 


that 
ees VE 
¢ Vii a (845) 
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This is the connection between the strength of the singularity at z=0 and the 
breadth 26 of the plate. We have shown now that linearized theory can be 
understood as a first order approximation of reentrant-jet theory. The same 
statement holds for Riabouchinsky flow (see [13]). 


9. Unsteady flow; method of the complex velocity 


The coordinate system and the units of length and time are chosen as in 
Section 8. It is assumed that the strength of the singularity at z=0 will vary 
harmonically with circular frequency w, and that the amplitude of the time- 
dependent part is small with respect to the time-independent part, but large 
with respect to second and higher powers of this latter. This last assumption 
enables us to treat the unsteady effects as linear perturbations of the linearized 
steady flow (second linearization). Hence all quantities can be split up addi- 
tively into steady and unsteady parts. Accordingly we put 


U=U+u with 4—%, ec, 


Vv=%U+y, De Oe, 
A=A-A,, has: A, dvs ah 
B=B,+ By, Need ae 
fe ale 
p=hyo+ hr b=pem. 


Remark. The time variable is denoted as usual by ¢ instead of by 7, differently 
from the notation used in Part I. Asa result of the unsteady motion the boundary 
condition on the cavity will change. For its derivation we start from the linearized 
Euler equations of motion for an ideal fluid: 


Ou CUns 4 ise? 
ra} ? 
ot Ox Q@ Ox (0.2) 
ov a ov 1 op 
ot Ox Q oy” 
On the cavity we have p=#,; thus - = 0, and by (9.2) we get 
Ou ou 
OE i we O (9.3) 
By use of (9.1) we obtain 
eae tam 
OU = 0. (9.4) 
The general solution of this differential equation is 
a, = ge 1°, (9.5) 


where g is an arbitrary constant, real with respect to 7, but possibly complex 
with respect to 7. 


338 J. A. GEURST: 


For w=w,+w, we now have the following mixed boundary-value problem: 
Determine w(z), holomorphic in the upper half of the z-plane, so as to satisfy 
the following boundary values on the real axis: 


= CHL HOO (0). Ol al tani 
020 21 wee ole. or shew ee (9.6) 
Pree ae, by, or) Im a==0. 


By means of the Schwarzian reflection principle this boundary-value problem 
can be transformed into a Hilbert problem with the following linear relations 
along the real axis: 
= 00 PO OO 
Oma aly wt + w= 2g", (9.7) 
[Ee GEICO LOE ig =O 
The solution is 


l are 
re —— a fp OX 
ee pene a fj Je (ht dx ig | pa B|/ 2 Sigs 
Z m yale “ Vx (I— x) (w—2) | a a Oe 


We now divide this expression for w into its steady part w, and its unsteady 


part w,. We find that 
| j2—l, Zi 
Wy = Ay — By ap. (9.9) 


Re 
NI pearreemecnaar a Med aden 
w= fe ‘Vale I) [5 
0 


x (ly—*) (¥—2) 


L 1 Ly z 1 ay a 
Ao 2 Vz(e—h) aes 2 \% 2 a | ent a eS ae 


For the relation between A, and the breadth 20 of the flat plate, the determination 
of By, and the connection between /, and the cavitation number o, we refer 
to Section 8. 

Assuming that in (9.10) the strength of the singularity at z=0 is given, we 
still need three other conditions for the determination of the four unknowns, 


(9.10) 


g, ie A, and Be These conditions read as follows: 


4. As in Section 8 we must have 


w=— pO alt Z=oo, 

or (9.11) 
w,=0 Au = Ca, 

This condition means that the perturbation velocities vanish at infinity. 


Phe LES Wi—= OP) Ate ze — 100% 
or, according to Section 8, 
res.W,=0 at =o. (9.12) 


This condition means that there is no source or sink at infinity. We already 
know that there is no vortex at infinity as a result of the symmetry of the problem. 
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As a result of this condition the pressure remains finite at infinity in the unsteady 
case too. 


3. Boundary condition (9.5) on the cavity has been obtained by considering 
the derivative of the pressure ~ along the boundary of the cavity. Hence it is 
possible that a time-dependent additive function, a constant being included as 
a special case, is not taken into account for the pressure in the cavity. In any 
case it is necessary to connect the pressure on the cavity with the pressure at 
infinity. This is the third condition. We now formulate it mathematically. 
According to (9.2) we have 


Cpa 0; OUo 
Ox age” 
and (9.13) 
Op ATR, ; Ou 
Bx ey e(fom+ ae 


The unsteady motion is not allowed to disturb the connection between the 
pressure at infinity and the pressure on the cavity, 7.e., we require that 


Pie fi ose OU: (9.14) 
or, by virtue of (9.13), that 
+0 
ae net)" = 
{ (Goat rh) ax =0. (9.15) 
Integrating, we obtain 
+0 
jo f adx=— %,(+0,0)=—, (9.16) 
or 
+0 2 
faaxa Ze. (9.17) 


Application of (9.11) and (9.12) gives 


Io 
“a —jox x a 
{3 a8 +B, =0, (9.18) 
# # (Iy—*) 
0 
IB . 

@ yma pn OE eda) 

T 2 Vx (j—%) Vx (Ip—*) 
i 2 9.19 
ee ee a hs 

Sey 1 Sie 0 af 2 t 


Since Ay,= By= — $29, it follows that 
Ts ark 285 
aS ee 
bf give ee dx +A,ly=0. 
Tt 
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(9.20) 
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By means of (9.10) and (9.20) we find that 


ly 
ey 1 pias VElo—8) ge oh fo . (9.21) 
ae n ein sink 8 (—x)} ()—2)* 


We can now apply (9.17) 


ly 5 
J oe ax 
—£ fe i ‘VE — é) ae [ V—#(—*) (§—*) 


0 Sg (9.22) 
iz pe dx (IS 
SO A oe et ee 


The inversion of the order of integration can be justified in the same way as 
in Section 5. We now use the following identities: 


é 
id? =e (4) ets arc cos (7 — ‘)), 


; (9.23) 
Hl aK Se 
L (-a(—*)t bo 
Integrating partially, from (9.22) we get 
Ip rs 
“a gal or - 
g dx oljo 
— ——— =Q. 9.24 
‘a Vx (lp—2x) a a 


By means of (6.16) the integral in this formula can be expressed in terms of a 
Bessel function : 
1 a0 ~ —1ig it \ 
<ohjo= pew ' Jo(w 2). (9.25) 
Remark. It is at once possible to write the solution (9.8) of the boundary- 
value problem (9.6) in the form 


, , 
Ree 1 e 1°* Vx(l—x) dx Ve , Zz 
ete Vz(z—D) i 4%—z eda z te | z—l~ (9.26) 


By putting A’= A, and B’= By, the conditions (9.11) and (9.12) are automatically 
satisfied. Finally it is possible to get from (9.22) to (9.24) without evaluating 
the first integral in (9.23). In this case one uses the fact that the integral can 
be considered as the solution of a certain boundary-value problem. 


10. Linearized unsteady cavity flow as a first order approximation 
We shall derive expressions for Q, and #, which can be compared in order 
to show that, for large w, #, can be identified with the first-order approximation 
to oe derived in Section 6f. For that purpose we express Q, in terms of g 
and w#, in terms of g. From (6.13) and (6.14) we find that 
ate e 1 ee = ete 
aaa Trechle Valea ge r Co 35 te a ‘fi 15 = ee 
0 


—o) 
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It has been shown that this expression is valid only for large w. Combination 
of (9.10), (9.20) and (9.24) gives us 


if LL oe ‘ 
o,= g 1 Ifesk aed dé 1, Jj 


m Vz(z—I,) ; 20 


; thy jviwe ae 
A(z—h)i J VE(Q—8) 


a |%> 


, (10.2) 


In Section 8 we have identified the following quantities, viz, ¢ and z, cy and Jy. 
If we may now identify g and ¢ also, the two expressions (10.1) and (10.2) become 
the same. If we put the strength of the singularity for #, at z=0 equal to the 
strength of the singularity for Q, at t=0, we infer from both formulas that 
indeed g and g must be identified. This gives the proof of the statement that 
for the unsteady case, too, linearized cavity theory is a first-order approximation 
of reentrant-jet theory for large w. From (6.22) and (9.25) we see also that ¢, and 
i, must be identified. 


Finally we derive an expression for the strength of the singularity at z=0 
in terms of # and w. This expression can most easily be derived from (6.13). 


According to that formula Q, has the following behavior near t=0: 


5 Dae - 
Q,~ —joh2*? Ya(— t)-3, (10.3) 
where (—?#)~# is positive for real ¢<0. Application of (6.6) gives us 
Agee a+2 2b ay 
Q,~—joh2t \2% a. (10.4) 


11. Unsteady flow; method of the complex acceleration-potential 


The coordinate system and the units of length and time are chosen as in 
Section 9. We start from the linearized Euler equations 


Ou Ou 1 Oop 

rie anrrmes 0 ox’ (14.4) 
ov ov 1 Op 

ot 0%, 6 OV 


In these equations we have linearized with respect to g,=1 the velocity on the 
cavity in the steady case. We introduce the function » by means of 


g=— PaPe : (44122) 
Q 
Hence (11.1) can be written as 
Ou Ou Op 


ot + Ox Ox 
ov ov O@ 


(F453) 


ot + Ox oy - 
@ is called the acceleration-potential. As a result of the incompressibility and 
irrotationality @ is a harmonic function. We can now introduce its conjugate 
harmonic function #, satisfying the relations 


EA ae 

Lovee (11.4) 
Cie So 

oy Ox 
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Now F(z)=@+7 is a holomorphic function of z=x-+7y for which we have 


oF aw éw (11.5) 


pe. Of tigers 
F(z) is the complex acceleration-potential. Now we may formulate the linearized 
cavity problem of Section 9 as the following mixed boundary-value problem 
for F(z). 
Determine F(z), holomorphic in the upper half of the z-plane and assuming 
on the real axis the following boundary values: 


—wo<x<0, y=0O or ImF=0, 
Oxia; g=0 or ReF=0O, (11.6) 
La2z<oo, y=0. or Imr=0. 


The general solution with integrable singularities is 


z—l Z 
Ez) a ae b les ; (11.7) 


where a and b are quantities dependent only on ¢ and real with respect to 7. 
Remark. This a is different from the a used in Part I. According to the 
second linearization we can write 


F=RK+R (with Rake, 
a= dg +a, a, =, ef”, (11.8) 


b= by +h, b= b,%. 


We identify Fj and wy. This is in accordance with the foregoing equations. Thus 
Z—l, Z 
ee ea ee 
For &, we now get the expression 


A 1 X 1 1 
De ee eee eae EY ee : 
: O72 AV Gr ° 2 Vz@—h) @— 
0 0 


e—l 
ba, |/ i De sas where A, =/1— 1. 


In this expression for /, three unknowns, v7z, 4,, a, and 0,, still occur. For their 
determination we have the following three conditions: 


Wa 


o 
4 


(11.10) 


1. F=0 at s=0, (41.41) 
v.e., the perturbation pressure vanishes at infinity. 

2. The velocity must be a single-valued function. Naturally this requirement 
is fulfilled for the horizontal component of the perturbation velocity u, by virtue 
of the symmetry of the problem; there is no vortex sheet in the wake. Thus 
Condition 2 expresses the fact that there is no line of discontinuity for the 
vertical component of the perturbation velocity v,. 

3. The unsteady part of the singularity for F(z) at z=0 and the given un- 
steady part of the singularity for w(z) there must satisfy an appropriate relation. | 


Unsteady Cavity Flows 343 


As to Condition 3, we shall prove that the unsteady parts of the singularities 
for F(z) and w(z) at z=0 are the same. 


According to (11.5) we have 


jo®, + — pce (11.12) 


from which it follows that 


Fp —jJwz 1) OF 
Wiel f eet Seat 
—00 (11.13) 


=F(2)—jo fd RO de. 


Use has been made of (11.11). Since the integral in the last member of (11.13) 
is finite at z=0, @, and & have the same behavior at z=0. q.e.d. 
Equating the strengths of the two singularities, we obtain 

olde haa etree j 
Gh sete = At 5, Ag a 
a) i} ° OVE 
This is the mathematical expression of Condition 3. Application of Condition 14 
gives us 


(41.14) 


4, tb; =0; (11.15) 


We shall now formulate Condition 2 mathematically. Since F(z) is a single- 
valued function, we can infer from (411.13) that the single-valuedness of #, is 
equivalent to 

f dF (0) dt =0, (11.16) 
C 


where C is a contour surrounding the cut from 0 to / in the z-plane. Since 
ay=Ay=— fa, we can write (11.10) as 


they meee Z—ly ; z O'My bacta! 
A) =a |/ z si loee pe Vz(z—1p) 


m ‘ (11.47) 
ol, z 1 _ hy 4 oi, d / z 
8 zZ—ly 2—ly Vz(z—1y) 4 d2 Vy 415 ~ 
where use has been made of (11.15). Condition (11.16) is now 
~ ess Ajo j f 2 
G1 i! é -dz+ SAI f gio 7 da=0, 14.18 
I) Ve@=I,) 4 = (eee 
C C 
or 
to Jox ava to i ie 
aly f SO ee Oe leah Lee): (11.19) 
Vx (ip—*) 4 ly—% 
0 0 
Replacing x by J,—% as variable of integration in the integrals, we get 
eis i lp Refol ie 
—jox a 
igh Ewa: ax ! oA] of ee |epeax=o, (11.20) 
F Vx (Ip—*) 4 F % 
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or, using (6.15) and (6.16), 


‘ Aj l I 
a, Jo(o 2) - oa [Jo [w 2 +4 Si (o2)| =0. (Ata) 
From (9.20) and (9.24) we can derive 
= | 
A, Jy (« 2) fae J( 2) +7f, (» 2)| == (0). (44-225 
Combination of (11.14), (41.21) y (11.22) yields 
eae and =k. (11.23) 


This means that both methods, viz, the method of the complex acceleration- 
potential and the method of the complex velocity, lead to the same results. 
This is important, because the method of the complex acceleration-potential is 
preferable in more complicated problems. 


12. The drag 


We calculate the drag by applying the momentum theorem to the interior 
of the contour indicated in Figure 5. C, represents the linearized version of the 


i SF Dicareenemmmernennemcentie le tel eile P52 RE 2 oy 
Us 
Gi 


Fig. 5. Contour for the calculation of the drag 


flat plate. The whole contour is denoted by C. We use the method of the 
complex velocity. 

It is clear by virtue of the symmetry with respect to the x-axis that the 
component of the force perpendicular to the main stream vanishes. Since the 
horizontal sections of C give no contribution to the x-component of the force, 
and since the contribution of the vertical end sections can be made arbitrarily 
small by making the contour thinner, the momentum theorem can be written as 


id = eh yeved 1+u)vdx, (12.4) 


where # is the pressure exerted by the fluid on the flat plate and the direction 
of C, and C is chosen positive, 7.e., such that the interior of C lies on the left. 


Since the drag D= a (p—p,)dy, we have 


eae 1+ u) ie ae i Strela Nes 


(12.2) 
--me[ S Sean ee 


By virtue of (8.15) 


Py ake) 
Wo u eR gh > at z=0, 
and by virtue of (40.4) 
B,~ —ijoh 222 me? —# at g=0 
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Hence we have 


Me ute e205, . U2 jot] 4 ie) 
w= Al saan [t+ jah 222 lz at z=0. (12.3) 


Substituting in (12.2), we find 


Bia 205s & ajo hE? dot), (12.4) 
and 
_— D 20 ‘ WR) aren 
Cp TOI Fae [A+ 2joh =? o |. (12.5) 


In (12.5) we have neglected the factor (1-++¢) in accordance with the order of 
approximation we have used throughout. In [13] we have given a formula from 
which it can be seen that for the steady case the insertion of the factor (1+) 
— already known as the (1-++a) rule — gives the correct next order of approxi- 
mation. 

13. Conclusions 


In the foregoing sections we have shown that a linearized cavity theory can 
be formulated in such a way that it becomes a first order approximation to the 
reentrant-jet theory. For the proof we have ‘chosen the simple example of a 
flat plate normal to the direction of the main stream since in that case the 
calculations are as simple as possible. The advantages of this choice are slightly 
diminished by the fact that we could compare our linearized theory with the 
first-order approximation to the unsteadily perturbed reentrant-jet flow only 
for rather large values of w, as discussed in Section 6. We had to make this 
restriction because the flat plate performs an oscillation in the direction of the 
main stream. For a lifting plate, oscillating normally to the main flow direction, 
the restriction concerning w drops out. With respect to the behavior of the 
flow at the end of the cavity, however, which is the important point in the 
steady case as well as in the unsteady case, we may safely assume that it is not 
influenced by the shape of the forebody in the special example chosen. 

We have assumed that the unsteady motion is small with respect to the 
steady part of the disturbance field caused by the body and the cavity. For 
the non-linear theory ,?.e., the reentrant-jet theory, this assumption was necessary 
in order to attack the problem in a reasonable way, and for the linearized theory 
it has enabled us to arrive at explicit results. Moreover we may expect that 
the condition of small unsteady perturbations will be satisfied in many operating 
conditions. The hydrodynamic stability of the cavity flow was not considered 
in this paper. It may be determined, however, from the results that we have 
obtained for the case of harmonic motion. 

We have shown that for the linearized cavity flow the method of the complex 
acceleration-potential leads to the same result as the method of the complex 
velocity, provided appropriate additional conditions are applied. In the case 
of the second method these conditions may be derived easily by comparing the 
linearized flow with the reentrant-jet flow. A condition which seems to be new 
is the one relating the pressure at infinity and the pressure on the cavity. From 
our analysis it is clear that in using the method of the complex acceleration- 
potential it must be required explicitly that the normal component of the per- 
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turbation velocity should be continuous in the wake. Since this condition has 
been overlooked hitherto, we mention it here. 

The linearization of the vertical flat plate into a singularity, which may be 
considered as a slight generalization of linearized theory as used thus far 
enables us to treat more complicated problems in a relatively simple way. E£.g., 
the effects of gravity on a cavity have been dealt with recently in that way. 
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Bounds on the Dissipation of Energ y 
in Steady Flow of a Viscous Incompressible Fluid 
around a Body Rotating within a Finite Region 


E. A. KEARSLEY 


Communicated by R. A. TOUPIN 


1. Introduction 


The problem of calculating the torque on bodies of revolution in a viscous 
incompressible fluid in steady flow has received considerable attention. Such 
calculations are needed in designing and calibrating rotational viscometers. With 
few exceptions, previous investigations have used STOKES’ equations of slow flow, 
in which the inertial effects of the fluid are ignored. The attempts made to take 
account of these inertial effects have usually depended on a perturbation technique 
in which deviations from Stokes flow are expanded in powers of some Reynolds 
number*. Such methods appear to be valid for very slow flows, but they suffer 
from the difficulty that the range of validity is not given by the method and 
can be determined empirically at best. 


When inertial effects can validly be ignored, so that SToKEs’ linearized 
theory applies, it 1s well known that dissipation of energy, and thus also the 
torque, can be bounded both above and below by certain functionals of a trial 
function **. What is more, if the true Stokes-flow velocity is selected as the 
trial function, these bounds converge, and the actual torque is determined. 


This paper explores the possibilities of extending these bounding techniques 
to the non-linear case where inertial terms are considered in full. The Stokes- 
flow solution is used as a trial function, and upper and lower bounds are found 
for the torque on a rotating figure of revolution in a finite fixed container in 
terms of the torque of the Stokes-flow solution, but only for the case that a 
quantity having the character of a Reynolds number is less than one. One 
valuable feature of the bounding techniques for Stokes-flow does not carry over 
in this generalization to the non-linear case; if the actual flow is used for a trial 
function in calculating the bounds given in this paper, the upper and lower 
bounds do not converge. The question of existence of a steady flow is not in- 
vestigated. 


* See for example D1 FrANciA: Bull. Unione Matematica Italiana, Ser. 3, 5 (1950). 
xk See Appendix I. 
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The principal results of this research are twofold: 

1. A rigorous maximum error of the Stokes-flow torque is displayed in terms 
of a uniquely defined quantity with the character of a Reynolds number. 

2. A technique is suggested for studying the hydrodynamics of highly non- 
linear steady flows. When equations are extremely difficult to solve, one is 
grateful for the information contained in inequalities. 


2. The Upper Bound on Dissipation 


The configuration to be studied consists of a figure of revolution rotating 
about its axis of symmetry and immersed in an incompressible fluid which is 
itself confined by a fixed outer wall; see 
Tig. 1. The region filled with fluid, assumed 
finite, will be denoted by R; the rigid rotating 
surface will be denoted by S,; while the fixed 
surface will be denoted S,. S will denote 
“\ So the boundary of R, the sum of S, and S,. 

The Navier-Stokes equation for an incom- 
pressible fluid can be written 


Piss = i(() 
(1) where 
Pig = — 0 Oi; + 2 d,;(U) — 04; %;. 
Fig. 1. A cross-sectional diagram of the In this and following equations, repeated 


geometry considered : s : ; : 2 
indices will indicate summation, and the 


usual tensor notation for Cartesian coordinate systems is used; # the pressure, 
ju the coefficient of viscosity, @ the density of the fluid, and 6;; the Kronecker 
delta. The 7‘* Cartesian component of the velocity vector wu of the fluid is 
written u;, and all subscripts following a comma indicate differentiation with 
respect to the corresponding Cartesian coordinate. In this notation, the rate of 
deformation tensor of the vector velocity field wu is defined by 


(2) d;;(U) = 3 (ui; +; ;)- 

By GREEN’s transformation, one can obtain 

(3) S 0; pi 48; =f Y;,; i, 4T 
s R 


where v is a vector field, continuous in R and continuously differentiable through- 
out the interiors of a finite number of volumes whose sum is R. The volume 


element in R is dt, and ds is the surface element directed toward the exterior 
Ot Ak. 


Equation (1) implies 
(4) Ae i348; = — J pv,ds;+ 2u fv; d,;;(u)ds,;— of v; u;Uu,;ds;. 
s s s 
But u-ds=0 on the impermeable bounding surfaces. We now specify that 


v=u on S so that equation (4) becomes 


(5) fv; Pi, 48; = 2u fu; d;;(u) ds;, 
Ss. S 
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and this quantity is just the rate of doing work on the fluid in R. We can call 
this quantity w M where w is the angular velocity of the surface S, and M is 
the torque acting on S,. 


Equations (3) and (5) with the symmetry of ¢ imply that 
(6) wo M= FM, 4 Pi =I di 5(0 v) p;, aT. 


By further restricting » to be solenoidal, i.e., Jv; ds;,=0 over any closed surface 
in R, we obtain 


(7) aM eye, OVE, ;j(u Ne as AO lnde U; d;;(v) at. 


Partial integration of the second integral on the right side of equation (7) 
yields 


(8) dee u, d; ;( CS NUS Be 
In virtue of the identity 
(9) Hhiepe U; i, a) HGH; Os =O, 


equations (7), (8) and (9) may be combined to give 
(10) OLE ARNE d;;(u Ue 20 Jt u,d;,(u) at. 


The restrictions placed on the vector field v up to this point are the following: 


4. Continuous in R and continuously differentiable in a finite number of 
volumes whose sum is F. 


2. Converges to u on S. 
3. Solenoidal in R. 


It is clear that we may choose v to be the solution to the Stokes-flow problem 
whenever the latter exists, although it is not necessary for w to satisfy any 
differential equation other than that its divergence shall vanish almost every- 
where, as is implied by its solenoidal character. 

It can be shown that 


(11) OM up aalM d;;(u) dt, 


which is just the statement that wM is equal to the dissipation of energy in the 
fluid. By analogy, we define My, the torque of the trial vector field, by 
(12) @ My = 2p Jas d;;(v) dt. 


Then, applying SCHWARz’s inequality to each term on the right of equation (10) 
with some manipulation results in 


(13) oM<o\/MM, +0¥|/7% 


20M 


fe WE 
J 
where ¥? is the maximum value of v; v; in R*. 


* When v is an axially symmetric Stokes-flow solution, it can be shown that it 
assumes the magnitude Yon the boundary of RF. 
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To find an upper bound for M, it remains to bound f wedt from above. 
To do this we consider GREEN’S transformation, K 


(14) fu; wi5 4s; Se eee Livin de: 
S 


where tp is any suitably smooth tensor defined in R and on the boundary Ss 
In particular, let y,;=4; x, +2; x,; then 


(15) f (tu; ; Xj + 20; 4; pH) ES — 4 J 0% x,d,,(u)dt+5 fwdr. 
S R 
The surface conditions on u may be used to reduce this to 
(16) f uty d55= 4 f sx d,(u)dt+5 fwdt, 
os R 


which, upon application of SCcHWARz’s inequality and some manipulation, yields 


(17) s[wrdr— 20 |/ futdr 72" s a 


where & is the maximum radius in R from the coordinate center, and we have 
used the definition 
(18) P= fue ds. 


Si 


Inequality (17) can be solved to give 


(19) | fears //22™ F(14 \) r 1+ 3 Soe). 


By combining this with inequality (13), one obtains 


a a 20% Sie ; 
(20) |Ms yM, +72” #yu(t | \/: te A 


It is easily seen that if 


(21) Nee 


40VR 


<7 

inequality (20) implies a finite upper bound on M. This upper bound can be 
calculated numerically for any particular geometry. With the loss of some 
sharpness, however, we can solve (20) explicitly. In Appendix II it is demon- 
strated that J is negative; therefore inequality (20) can be shortened to give 


(22) Ms ee when Np<1. 


By observing that (1++-a)'<1 os one can sharpen this result somewhat to give 


(23) yu < VM 4 _5utNe 

1—Ne ' 302M, 
Undoubtedly even sharper approximations to the solution of inequality (20) 
can be obtained, but further efforts in this direction seem hardly worthwhile 
in light of the fact that it can be solved numerically for any particular geometry. 
Inequalities (22) and (23) are rigorous upper bounds on y™. 
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3. The Lower Bound on Dissipation 


A lower bound may be obtained more briefly. We choose a trial tensor , 
symmetric and such that 


(24) Yii,7=%i mR, 


where g is a bounded scalar field in R. Tensors satisfying this requirement can 
be simply constructed in the following way*. Choose a symmetric tensor field, 
say a, with finite second derivatives defined in R. Then 


(25) Vii = Eikl Ejmn Linnm + 4 0%; 


is a solution of equation (24), where ¢;,,, is the usual permutation tensor density. 
By GREEN’s transformation, one can write 


(26) Suzy; 4s; = fu; piydt + fugit. 
S R R 


Because of the symmetry of and the boundary condition on wu, equation (26) 
can be written 


(27) J) U5 Piz TS; al d;;(U) yij aT. 
Further, by SCHWARz’s inequality, we have 

LUG Yi j dsj) 
(28) oM=2u { d;;(u) djj(u) de 2 2n— 


Inequality (28) constitutes a calculable lower bound on the energy dissipation 
in R. In particular, we may choose for the trial tensor \ the rate of deformation 
tensor of the Stokes-flow velocity field satisfying the prescribed boundary con- 
ditions, v2z, 

Se seed 

(29) Wij =4,;(V), 9 2m 


where ?’ is the pressure field associated with the Stokes-flow solution. Then 


[ fv dis (v) asi) 


> i | 
60) OM SAH Fase) dae) tt 
and hence 
(31) M2 My). 


4. Conclusions 
1. It is interesting to observe that the greatest lower bound contained in 
inequality (30) is inequality (31). This fact is clearly implied by inequality (41) 
of Appendix I. In general, this analysis is too rough to expect convergence of 
the upper and lower bounds for any trial functions, a fact which is in contrast 


* See Dorn & SCHILD: Quart. Appl. Math. 14, 209—213 (1956) and TRUESDELL: 
Arch. rational Mech. Anal. 4, 1—29 (1959). 
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to the Stokes-flow case discussed in Appendix I. An interesting question is 
whether a modification of this bounding technique could be sufficient to assure 
convergence of the bounds. 

2. The upper bound on the dissipation is lost if Np is greater than one. In 
a later paper, it will be shown that this fact is a result of the roughness of the 
approximations employed in the analysis and has no hydrodynamic significance. 

3. Inequalities (143) and (17) are enough to specify an upper bound on total 
kinetic energy of the flow. A simple and less sharp form of this bound is 


8022 a0 M, 

2 < 

(32) g fu dts 254. “(t= Ng)? © 
R 


which is valid when Np is less than one. 


4. The quantity Np is a dimensionless number which may be considered a 
Reynolds number, uniquely defined and calculable for all cylindrically symmetric 
cases and all cases where the Stokes-flow solution is known. It serves as a rigorous 
criterion for the error introduced in calculating dissipation by ignoring the inertial 
terms in the Navier-Stokes equation. Put another way, it specifies conditions 
which guarantee that the dissipation calculated from a Stokes-flow solution will 
be within a specified range of the dissipation according to the full Navier-Stokes 
steady flow solution. 


5. Although this analysis has been done for the case where the inner surface 
rotates and the outer surface is fixed, it is clear that a simple extension of the 
analysis can take into account rotation of the outer surface as well. 


Appendix I 


The Stokes-flow equation for an incompressible viscous fluid can be written 
(33) 937, 7=0, where Vii = — £0;;+ 2ud;;(u). 


If uw is the solution to this equation and v is a smooth divergenceless vector 
field* equal to u on S, then by GREEN’s transformation 


(34) JS 0; 9:5 48; = J 0; 9:47, 
Ss R 


where © and w satisfy the boundary conditions of the previous problem. By 
equation (33) this becomes 


(35) wo My= 2 f dj;(u) d,;(v) dr. 
R 
An application of ScHwarz’s inequality to equation (35) yields 


(36) (o My)? S 4u2 f d;;(u) d;;(u) dt f d;;(v) d;;(v) dt. 
R R 


* More precisely, » must satisfy the three conditions written out previously 
following equation (10), 
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Finally, in light of equation (5), we have 
(37) o My S 2u f d;;(v) d;;(v) dt, 
R 


which constitutes an upper bound* on w My. 


For a lower bound on the dissipation of energy, we select as a trial function 
a smooth symmetric tensor field, , which is a solution of 


(38) Dig = | 45 


where g is a continuously differentiable scalar field in R, converging to a finite 
value on S. Then by GREEN’s transformation 


(39) J up) 4s;= Su, jppidt+ Jug at. 
s R R 
But in view of the boundary condition on wu and the symmetry of , 
(40) fin d= 8S 45(4) ij ae. 
A R 
Then by SCHWARz’s inequality and equation (5) we have 


( fu; viz dsj), 
41 (OA es 
a Neve Ra Yan ee 


This inequality should be compared with inequality (28). 


Inequalities (37) and (41) constitute respectively upper and lower bounds of 
the dissipation of energy associated with a Stokes-flow solution of the problem. 
In particular, if the trial vector v and the trial tensor are taken as the Stokes- 
flow velocity and rate of deformation tensor, respectively, these bounds will 
converge to the dissipation associated with wu. 

Although these bounding principles for the linear Stokes equation seem to 
be well known, I have been unable to find in the literature a good reference for 
the upper bound. Almost identical, however, are the principles of DrRICHLET 
and THomson which bound the Dirichlet integral of solutions to the Laplace 
equation **. 


Appendix II 
The integral J has been defined in equation (18), vzz: 


(18) Te aed sa 


Sy 


* This bound can be recognized as a form of the well known Helmholtz Theorem; 
see for example SERRIN: Handbuch der Physik VIII/1, p. 258. Berlin-Géttingen- 
Heidelberg: Springer 1959. 

*xx See for example J. B. Diaz: Upper and Lower Bounds for Quadratic Integrals, 
and at a Point, for Solutions of Linear Boundary Value Problems. Proceedings of 
a Symposium on Boundary Problems in Differential Equations, Madison Wisconsin, 


1959. 


; te gl a", +5 et 7, f' Wise ee 
ied cake 1 al rare 
ah are 

; 


354 E. A. KEarsLtey: Torque Exerted by a Viscous Fluid . 


The velocity on the surface S, is known, since S, is the surface of a rigid 
rotating body with angular velocity specified. If w is the angular velocity vector 
of S,, we have 
(42) U; = E540; %, On Sy. 


Let us consider equation (42) as extending the vector field uw into the interior 
of the volume bounded by S,. Then by GREEN’s transformation, noting that w 
is a fixed vector, we get 


(43) T= —5 f &:54 &j1m Oj Oy Xp Xp AT, 
e 


where R’ is the region bounded by S, and the negative sign is due to the fact 
that the outward drawn surface element of the region R is the negative of the 
outward drawn surface element of the region R’. We observe, finally, that 


(44) I=—5fwdrs<o, 
= 
which proves I to be negative. 
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